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ABSTRACT 


In this report an analytical method is presented for predicting 
lateral-directional aerodynamic characteristics of light twin-engine 
propeller-driven airplanes. 

This method is applied to the Advanced Technology Light Twin- 
Engine (ATLIT) airplane. The calculated characteristics are correlated 
against full-scale wind tunnel data. 

The method predicts the sideslip derivatives fairly well, although 
angle of attack variations are not well predicted. Spoiler performance 
was predicted somewhat high but was still reasonable. The rudder 
derivatives were not well predicted, in particular the effect of 
angle of attack. The predicted dynamic derivatives could not be 
correlated due to lack of experimental data. 
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, -1 .-1 

deg , rad 


, -1 ,-l 

deg , rad 


■ -1 ,~1 
deg , rad 


, - 1 A- 1 

deg , rad 


(AC Y 3 ) n (Ai) 


(AC ) 

n( V 


Increment in nacelle contribution to C„ 
due to increased dynamic pressure 
Increment in nacelle contribution to C^, 
due to propeller sidewash 


*3 


. -1 .-1 
deg , rad 


,, -1 

deg , rad 



! 1 
I r„ 


Symbol 


“S’, 


power 


LIST OF SYMBOLS (continued) 
Definition 

Increment in Cy due to power 

8 


Dimension 


Dynamic pressure increase behind the 
propeller as a ratio of the free stream 
dynamic pressure 


Dynamic pressure increase at the nacelles 
due to propeller as a ratio of the free 
stream dynamic pressure 

Spoiler lift effectiveness in terms of 
the change in zero-lift angle of attack 

Rudder deflection 
Spoiler deflection 


deg, rad 


deg, rad 


deg, rad 


Variation of viscous drag with angle 
of attack 


Downwash gradient behind the propeller 


Upwash gradient at the propeller 


Variation of sidewash with wing tip • 
helix angle 


Variation of sidewash with sideslip 


Variation of propeller sidewash-with 
sideslip 

Downwash angle behind the prooeller 


Son-dimensional spanwise station, 
y/ib/2) 


xx ii 
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LIST OF SYMBOLS (continued) 


1 

Symbol 

Definition 

Dimension 

Tt 

e 

1) Wing twist 

deg 

ft 

d 


2) Effective spoiler sweep angle 

deg 


A 

Sweep angle 

deg 

Si 

n 

A s 

Corrected sweep angle, 
tan -1 (tan(A c ^)/B] 

deg 

Lj 

X 

Taper ratio 


M 

y 

^TE 

Section trailing-edge angle 

deg 


Subscripts 




c/2 

Half-chord line 



c/4 

Quarter-chord line 



c.g. 

Center of gravity 



e 

Exposed panel 


'.} 

f 

Fuselage 



fh 

Fuselage-horizontal tail combination 



h 

Horizontal tail 



/ e 

Leading edge 


:] 

M 

High subsonic nach number 



max 

Maximum 



n 

Nacelle 



prop off 

Propellers removed 



power on 

Propellers installed and power on 



r 

Root 



t 

Tip 



V 

Vertical tail 





LIST OF SYMBOLS (continued) 


Symbol 

Definition 

w 

Wing 

wf 

Wing-fuselage combination 

wfh 

Wing-fuselage-horizontal tail 
combination 


xx iv 
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Tills report describes work carried out under the second phase 
of a project performed by the Flight Research laboratory of the 
University of Kansas sponsored by Grant N5G 1574 from the National 
Aeronautics and Space Administration (NASA), Langley Research Center. 
The purpose of this project was to correlate theoretically predicted 
aerodynamic characteristics of the Advanced Technology Light Twin 
(ATLXT) airplane with full-scale wind tunnel and flight test data. 

The original phases of this project were 

1. correlate theoretical predictions of longitudinal 
aerodynamic characteristics with full-scale wind 
tunnel data; 

2. correlate theoretical predictions of lateral-directional 
aerodynamic characteristics with full-scale wind 
tunnel data; 

3. correlate the results of point 1 and point 2 with 
flight test data. 

The results of the longitudinal correlation (point- 1) are 
presented in Reference 1. This report deals only with the lateral- 
directional correlation (point 2). Full scale wind tunnel data 
used in the correlation were taken from Reference 2. The methods 


used in the prediction were taken from References 3 and 4. 

The purpose or this report is to present an analytical method 
for predicting the lateral-directional aerodynamic characteristics 




1 



full-scale wind tunnel tests to identify areas where current 
predictions fall short. 

This report was originally to correlate the predicted char- 
acteristics with flight test data. However, aside from some 
performance data presented in Reference 5 , sufficient flight test 
data i3 not yet available to perform a meaningful correlation. 
After consulting with NASA, it was agreed to delete the f.’.ight 
test correlation phase. Additional work of a similar nature 
(full-3C3le wind tunnel correlation) is presently being performed 
on another airplane in place of the ATLIT flight test correlation. 
Some derivatives for which no wind tunnel data were available 
are still presented for future comparison when this data becomes 
available. 
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CHAPTER 2 


THE ATLIT AIRPLANE 

The ATLIT airplane is a Piper PA-34-200 Seneca I general aviation 
low-wing monoplane wich a number of high technology modifications. The 
major physical characteristics are listed in Table 2.1, and a three- 
view drawing is shown in Figure 2.1. 

The improvements implemented on the ATLIT were 

1. The original untapered, aspect ratio 7.25 wing having 
a 65,412 airfoil and an area of 19.4 m' (208.7 ft') 
was replaced by a new tapered wing (taper ratio 0.5) 
with a higher aspect ratio (10.32) and a smaller 

area (14.4 m-). The new wing used a 17-percent thick • 
GA(W)-1 airfoil. 

2. Full span 30-percent chord Fowler flaps replaced 
the original 20-percent chord plain flaps. 

3. Spoilers were in.-: railed for roll control instead of 
ailerons . 

The new wing area was made smaller and the aspect ratio increased 
to provide better drag characteristics Chan the original wing. The 
GA(W)-1 airfoil was used for its high lif t-to-'drag ratio and high 
maximum lift coefficient. For acceptable landing speeds, full 
span flaps were used with spoilers providing roll control. 

During flight test, supercritical propellers were also evaluated. 
However ,• these were not used in the full-s'cale wind tunnel tests 
and- are therefore not considered here. 


In addition a number of changes were made 
tunnel test to clean up the drag of the ATLIT. 

as built' and "fully clean" wind tunnel tests, 
in this report are done using the "fully clean" 


'n part of the wind 
This resulted in 
All comparisons 
wind tunnel data. 
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Table 2-1: Specifications 


of the ATLIT Mr plane 


Winn ; 

Are*. Q : ( {c -) 

Spaa, a (( C ) 

Aapect redo 
Thiclxesa ratio 
Dihedral, dag 
Tapac ratio 

Incidence angle at root. J (g 
Incidence aogl. at dp, deg 
Leading «d g « , nd trailing ad ?c 
Mean aerodynamic chord, a (( z) 
Aoot chord, a (In) 

Tip chord, a (t n ) 

Airfoil 

Fljp : 

Type 

Area (total). a : (f.2j 
Span/aide, a ((.) 

Chard, pert, of wing chord 
Inboard wing station, a (in) 
Outboard wing station, n (in) 
Leading edge retracted, pore, of 
Maximum deflection, deg 
Maximum travel. „ ( ia ) 

Soot chord, m i in) 

Tip chord, m ;in) 


14. AO (155.0) 

12.19 (40.0) 

10.32 

0.17 

7 

0.5 

0.5 

-2.5 

»««P. deg 3.67 

1.225 (4.013) 
1.575 (62.0) 

0. *87 (31.0) 
CAOO-1 

Foul.c 

' . > 3.56 (33.3) 

5.15 (16.91) 

30 

0-711 (23.0) 

5-567 (231.0) 

"ing chord 70 

40 

0-343 (13. 50) 
0-445 (17. iO) 

744 (9-n2> 


0 



Table 2.1 • 


Continued 


Sootie r- 

Type 


Ar * 4 (total), a 2 (f c 2 , 

Span/alde, a ( Ia) 

Inboard wing station. a ^ 
Outboard win, .ration, Q (|o) 

Hing *’ ?erc - «f Wn* chord 
!<axi=ua deflection, deg 
Kooc chord, ei ( In) 

Tlp a (in) 


Tail ; 

Stabiiacor : 


Arej (Including tab) , {fc .’j 
s P«n. a (f C ) 

c "°rd (constant), a {ln) 

Aspect racio 


Suce P angle, deg 


Dihedral, deg 

Hln ? e line, pert. 
Air:‘oj,x 

Vertical; 


° c * chord 


Area (including tab), „ 2 ^. 2 , 
s Pan, a (fz) 

Asoec t ratio 1 
Tap-.-r racij 
R,, ot chord, a (l„) 


Diffarential'. 

Upper Surfaci 

0.488 (3.2J, 

1.21S (127.0, 

■ 2-365 (101.0) 
3 -2?l (1:3.0, 
70 
60 

°-I77 (6.973) 
0,1 24 (4.880) 

(38.7) 

4 -13 (13.36) 
'0021 04.29) 
4.73 
0 
0 

26.6 

N'ACA 0010 

*•35 (19.9) 

1.S2 (3.00) 

1.19 

0.42 

l "'03 (67. 33, 


■Slat- 


o 



Table 2.1: Continued 


Kean aerodynamic chord, a (in) 

Leading edge sweep angle, deg 
Airfoil 

Stabilator cab: 

Area, a 2 (ft”) 

Span, a (in) 

Chord (constant), a (in) 

Tab hinge line to stabilator hinge line, 
a (in) 

Rudder: 

Area (including tab) , a 2 (ft”) 

Span, a (in) 

Chord (constant), a (in) 

Rudder tria tab: 

Area, a 2 (ft”) 

Span, a (in) 

Chord (constant), a (in) 


1.231 (30.47) 
39.92 
NACA 0009 

0.539 (5.3) 
3.023 (119.0) 
0.173 (7.0) 

0.462 (18.17) 

0.706 (7.6) 
1.525 (60.02) 
0.449 (17.69) 

0.153 (1.7) 
O.S95 (35 25) 
0.173 (6.9) 


Fewer? 3 ant : 

Number or engines 

Manufacturer 

Model 

Takeoff rpa 
Takeoff power, hp 


2 

Lycoaing 

IO-360-CIE6 

2700 

200 


Propeller : 

Manufacturer 

Model 

Number of blades 
Otamecer , ra (In) 


Hacrrell 

HC-C2YK-lCZF/:< 


1.531 ( 7 h.O) 


7666A 



Table 2.1 


Concluded 



-1 

i 

■j 

' 3 

i 


\ 



V! 

1 


R 

& 


O 


Hi 

i-i 

u 


r-> 


a 


\A 

& 


i-i 

> r 

y 




■ i 
0 



» 
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'I 


•1 


i 


i 


t 
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Ualghca and loading : 

Cvaaa uaighc, H (lb) 

Etspcy weigh:, H (lb) 

Uaelul load, N (lb) 

Ulng loading (at groaa velghc), 
X/a 2 (Ib/f t*) 


17,792 (4,000) 
11.01S (2,477) 
6,774 (1,523) 

1,236 (23.8) 
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r? 








2.1 Geometric Parameters of Wing and Tail 


Depending on which characteristic has to be determined, the total 
area or the exposed area of the wing and the tail are considered. The 
total planform is considered to extend through the nacelle and the 
fuselage, while the exposed planform terminates at the fuselage. 
Pertinent dimensions for the wing, the horizontal tai! and the vertical 
tail are shown in Figures 2.1.1, 2.1.2 and 2.1. j, respectively. 

Table 2.1.1 lists the geometric parameters of the wing and the tail 
pertinent in the analysis. 
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Horizontal Tall | 

Vertical Tell 

Total 

Exposed 

£x(>osc<i 

3.60 (33.7) 

3.25 (34.9) 

1.75 (10.8) 

4.13 (13.56) 

3.73 (12.23) 

1.52 (5.0) 

4.75 

4.28 

1.33 

0.871 (34.29) 

0.871 (34.29) 

0.723 (28.45) 

0.871 (34.29) 

0.B71 (34.29) 

1.575 (62.0) 

1 

9 

0.459 

0.871 (34.29) 

0.871 (34.29) 

1.201 (47.30) 

1.033 (40.68) 

0.932 (36.69) 

0.334 (13.15) 

0 

0 

- 

0 

0 

40.0 

0 

0 

! 34.5 

0 

0 

29.0 






























Figure 2.1. A 


Dofi.ni.Cion sketch of the fuselage and nacelles 



CHAPTER 3 




PRESENTATION OF RESULTS 


In this chapter./the predicted lateral-directional characteristics 
of the ATLIT airplane will be compared with full-scale wind tunnel 
data of Reference 2. 


3.1 Sideslip Derivatives. C„ , C , C. 

— * Yp.- 1 — 


The methods used in predicting the sideslip derivatives of the 

ATLIT airplane are presented in Sections 4.1 to 4.3. In Figures 

3.1.1 to 3.1.6 the predicted derivatives C , C , and C, are pre- 

3 n 2 3 

sented and compared with full-scale wind tunnel data for boch propellers 
off and power-on. 

Figures 3.1.1 and 3.1.2 show the propeller off and power-on 

Cy . In general the results show good agreement, althouen the 

3 


effect of angle of attack has not been properly accounted for — 
especially for the power-off case. It is felt that the angle of 
attack effect on the fuseJage and fuselage-vertical tail combination 
has not been accounted for. 

The C of the full-scale wind tunnel test shows a strong 
n 3 

variation with angle of attack which the theoretical methods did 

not predict. While the average values agree fairly well, the 

predictions show a basically linear change >n C with angle of 

n 3 

attack. With power on, the prediction is even worse, as the wind 

tunnel C varied even more with angle of attack. As is discussed 
r '3 

in Section 4.2.2, the wing position plays an important role in the 

fuselage contribution to C . However, there is not enough' data 

n 3 

available to account accurately for this factor. 
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Figures 3.1.5 and 3.1.6 show the propellers off' and power on 

C. prediction compared with wind tunnel data. For the propellers 
8 

off case the prediction is somewhat high, but the variation with 
angJ.e of attack agrees with wind tunnel results. With power on, 
the wind tunnel tests show a more pronounced variation with angle 
of att<- .L At the low angles of attack the prediction is not 
very good, but it improves as the angle of attack is increased. 

3.2 Control Derivatives 

The control derivatives for the spoiler and the rudder are 
predicted in Sections 4.4 and 4.5. Figure 3.2.1 presents the 
spoiler effectiveness for a = 0 as a function of spoiler deflection 
The predicted effectiveness is slightly higher than the wind tunnel 
data. Nevertheless the prediction is considered quite good, con- 
sidering the general nature of the prediction method and the lack 
of data for Chis class of airplane. 

Figure 3.2.2 presents the yawing moment due to spoiler deflec- 
tion. Again the prediction is slightly higher than the 'wind tunnel 
data, although the overall trend is predicted quite well. 

In Section 6.1.5 it is shown that the effect of power on 
the spoiler characteristics were significant in the wind tunnel 
tests; no method was found to predict this effect. 

Figures 3.2.3 through 3.2.5 compare the predicted rudder 
derivatives with full-scale wind tunnel data. In general the 
predictions are net very good. The predictions show little 
variation vita angle of attack ('except, of course, for C. ), 
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while the wind tunnel data exhibits considerable variation with angle 
of attack. There is little variation in the rudder derivatives with 
power in the wind tunnel data. It was assuned that there was no 
power effect on the predicted rudder derivatives. 

3.3 D ynamic Derivatives 

The dynamic derivatives of the ATI. IT airplane were originally 

predicted for comparison with flight test aat--,. This data is not 

yet available; however, the data is presented here for completeness. 

Power effects were only considered significant for . figures 

P 

3.3.1 through 3.3.5 show the predicted derivatives as a function 
of angle of attack. 





'• 'Pi'urt'os ■■ 

terar ison of predicted C,. with full scale wind 

i , 

tunnel data (propellers removed , N, = 


1.3 Million) 

















a Degrees 


Figure 3.1.5: Comparison of predicted C. with 

a ■ 

full scale wind tunnel data (propellers 
removed, N_ = 2.3 million) 


3 






















WIND TUNNEL 
PREDICTED 



a 'v Degrees 


Figure 3.2.4: 


Comparison of predicted C with wind tunnel data 









n 




4 


6 


8 







guro 3.3.J: Prediered C or the ATI. IT 
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CHAPTER 4 
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PREDICTION OF PROPELLER-OFF STATIC STABILITY 
AND CONTROL CHARACTERISTICS 

In this chapter the propeller-off static stability and control 
characteristics will be discussed. The methods of Reference 3 are 
used for most all predictions. However, since Reference 3 does not 
include methods for predicting spoiler control characteristics, 
Reference 4 will be used. 

4 . 1 Side Force Derivac ' -e , C 

d 

The side force due to sideslip derivative, C y , of the Complete 

S 

airplane is found by considering the contributions of the following 
components : 

(1) Wing, including dihedral 

(2) Fuselage, including wing-fuselage interference 

(3) Nacelles 

(4) Vertical tail, including the interference of the wing, 
fuselage, and horizontal tail. 


4.1.1 Wing >.ntribution, (C. 


yj + (c v,>r 
J r=o 


The wing contribution to C,. with no dihedral is found from 

6 

Equation 4. 1.1.1 from Reference 3. This expression is valid for low 
subsonic mach numbers and includes coinoressibilitv effects. 


(C., ) 


6 tan (A .,) sin (A 
c/4 w 


c/4 'v 


?/w r = C 57.3 ta [A 3 n 
if w v 2 


;os (A ] 

C/ v 


(4. 1.1.1) 
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a 


i 




u 


i 





si 






a 


where 

C T is che wing alone lift coefficient obtained from Figure 4. 1.1.1 
L 

W 

is the aspect ratio of the wing 

is the quarter-chord sweep of the wing 

B, “ Vl-M 2 cos 2 A ,, 

2 c/4 

where M is the Mach Number. 

For typical general aviation aircraft with essentially straight 

wings and noderate aspect ratio (A > 6) , this contribution is insig- 

w 

nificant. For the ATLIT airplane 

(C Y ) ■ 2.32 x 10 -7 C 2 per deg. 

X B r w 

The increment in side-force derivative due to wing dihedral nay be 
approximated at low subsonic speeds by the following equation from 
Reference 3s 

(C )„ =■ -0.0001 T <4. 1.1. 3) 

For the ATLIT airplane che wing dihedral angle 7 is 7.0“. Therefore, 

(C Y ). = -0.0007 per deg.' 

3 1 • 

The summary calculation of the wing contributions to C,, are given 

3 

in Table 4. 1.1.1. 

4.1.2 "Fuselage Contribution to C y 

3 . 

The fuselage side-force due to sideslip contribution is the result 
of the side forces produced by the body and the wing-bodv interference 
effects. The fuselage alone is the main contributor ana is affected by 
the size and shape of the body. The wing-body interference contribution 
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is considered Co be only a function of wing location. The total 

fuselage contribution to C„ at subsonic Mach numbers is given by 

8 

Equation 4. 1.2.1 from Reference 3 (based on S ) . 

_2/3 

( C Y ) f ■ ) 2/3 — * per deg (4.1.2. 1) 

8 3 f— w 

v 

where 

(C Y ) is the contribution of the fuselage only based on 

8 ' 
v 

two-thirds of the fuselage volume and is equal to but opposite in 
sign to the potential flow part of the fuselage lift curve slope from 
Section 4.3 of Reference 1. 

is the wing-body interference factor obtained from Figure 4. 1.2.1 


Table 4. 1.2.1 summarizes the effects of the fuselage on C y . For 

V 

the ATLIT airplane 

(Cy ) f = -.00299 per deg. 

8 

4.1.3 Nacelle Contribution to C,, 

* ^ 


Due to a lack of definitive methods for calculating the contri- 
bution of the nacelles to C ‘ , Reference 3 presents Che following 

3 

empirical procedures: 

(1) The effective nacelle length is considered to extend only 
to the wing leading edge (see Figure 2.1.4). 

(2) The nacelle contribution is approximated from Equation 4. 1.3.1 
for bodies of circular cross section (based on 5 ). • 

M 


(C v )• 


2(k 0 -k. ) (S ) 

2 lx max 
n 


57.3 S 


per deg 


(4. 1.3.1) 



where 


n is the number of nacelles 
n 

is the effective cross section area equal to a circular cross 
' n 

section with a diameter equal to the maximum depth of the nacelle, d R , 
from Figure 2.1.4. 

t * ie re< ^ UCCl ^ niass factor obtained from Figure 4. 1.3.1 
as a function of effective length and maximum depth. 


(3) To account for flow interference effects, the (C^ ) n predicted 

3 

by Equation 4. 1.3.1 is reduced by one-third. 


(C Y ) = -2/3 n 

Yg n n 


2(k 2 -h 1 )(S ;c ) mgx 

n 

57.3 S 


per deg 


(4. 1.3. 2) 


The calculations for the ATLIT airplane are summarized in Table 

4. 1.3.1. The total nacelle contribution is 

(C Y ) = -.00048 per deg 

x 3 n 


4.1.4 Vertical-Tail Contributions to C,, 


The contribution of ; .e vertical tail to Cy is affected by the 

3 

location of the horizontal tail, the fuselage crossflow at the tail, 
and the wing-bodv induced sidewash. 

The horizontal tail can serve to increase the loading on the 
vertical tail by the so-called "end plate effect" when the horizontal 
tail is at a relatively high or low position. For mid span positions 
the effect of the horizontal tail is insignificant. 

The effect of the fuselage in sideslip is to increase the lo'-ai 
cross-flow velocity across the top of the body. This increased velocity 
causes an increase in vertical tail effectiveness. 


38 



f? 


; As Che sidewash due Co a yawed wing is small, Che primary sidewash 
coneribucicn is from che body and arises from che body vortex system 
produced in sideslip. This vortex system produces lateral velocity 
components which affect the vertical tail. 

The method of Reference 3 accounts for the horizontal tail and 
cross flow factors by computing an effective aspect ratio. The effective 
aspect ratio of the vertical tail is 


eff 


a&2L) 

v A 


1 + K h t-x;-- » 

v(f) 


(4. 1.4.1) 


where 


A y is the actual geometric aspect ratio of the vertical tail 


V(f) 


li 


u 


1 


% 


rt 

• : l 


is the ratio of the aspect ratio of the vertical tail in 


the presence of the body to that of the isolated vertical tail, ob- 
tained from Figure 4.1.4.1(a) 


v(fh) 

A v(f) 


is the ratio of the aspect ratio of the vertical tail in 


r-5 


the presence of the horizontal tail and the fuselage to the aspect 
ratio of the vertical tail in the presence of the fuselage 'only , ob- 
tained from Figure 4.1.4.1(b) 

is a factor which accounts for the relative size of the 
vertical and horizontal tails , obtained from Figure 4. 1.4. 2 

Table 4.1.4.1(a) shows the summary calculations to obtain the 
effective aspect ratio of the ATLIT airplane. It shows an effective 
aspect ratio of 2.46 compared to a geometric aspect ratio of 1.19. 

The lift curve slope of the vertical tail is found using the 
standard Polhamus equation (4. 1.4. 2) by using the effective aspect 
ratio (based on effective vertical tail area, S ) 
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n 




I 


xs . 


1 


3 


TT A 


(C L ^v(hf) 
a 


eff 


(4. 1.4. 2) 


2 + 


Vfr IB 1 J+tan2(A «« , '' ,+1 


k 2 

V 


where 


< C / >v 
a 

2ir 


(4. 1.4. 3) 


(C. ) is the section lift curve slope of the vertical tail obtained 
‘a 


from Section 4.1 of Reference 1. 


Bj 2 » 1 - M 2 


M is the Mach Number 


(4. 1.4. 4) 


(^c/2^v * s m ^-chord sweep of the vertical tail. 

Table 4.1.4.1(b) shows the summary calculations. For the ATLIT 
airplane (based on S v ) 

(C L ^v(fh) = ‘° 494 ?er deg 
a K k ' 


:v>\j 




The total vertical tail contribution to C y is given in Eauation 

3 

4. 1.4. 5 from Reference 3. It modifies the vertical tail lift curve 
slope to account for wing wake and body sidewash. 


(C Y J v(wfh) = " k l (C L ) v(fh) (1 + 36 ) - S 
3 a q w 

* ‘OB 


(4. 1.4. 5) 


where ' 


■ is. factor which accounts for the relative size of the body near 
the vertical tail to the size of the tail, from Figure 4. 1.4. 3. 


, .q, S /S 

(1+|~) — = 0.725 + 3.06 — 7 ,7 

33 - x + cos (..c/4) 

q 'v 


0.4 


7— :— + 0.009 A 
(w.) w 

t-j 


(4. 1.4. 6) 


.9 

i£l ‘ 




(A c /^) y is the quarter-chord sweep of the vertical tail 
z is the vertical distance from the centerline of the equivalent 

W 

fuselage to the quarter-chord point of the root chord of the exposed 
wing panel, obtained fron Figure 2.1.4. 

(w^) y is the depth of the equivalent circular fuselage at the wing, 
obtained from Figure 2.1.4. 

Table 4.1.4.1(c) shows the summary calculations used to find the 

vertical tail contribution to L . For the ATLIT airplane (based on S j 

6 W 

(C Y 3 ) v(wfh) " - 0 ’ 0056 per de 3 


4.1.5 C of the ATLIT Airplane 
3 

The total side-force derivative of the complete airplane (propellers 
off) is 

\ ' “V^-O + (< V r + “V* * “V" + <A - 1 - 5 - 1> 


Table 4. 1.5.1 summarizes the contribution of each component. For the 
ATLIT airplane 

(C„ ) =* -0.0098 per deg 

prop 

3 off 

The predicted and experimental C„ is shown in Figure 4. 1.5.1. In 

3 

general it shows fairly good agreement, although it is obvious that 
the effect of angle of attack has not been properly accounted for. 

These effects are probably in the wing-body interference and sidewas'n 
contributions. 



Symbol 

Description 

Reference 

Magnitude 

K 

Mach numb a c 

- 

0.081 

A 

V 

Wing aspect ratio 

Table 2.1.1 

10.32 

(A c /Jv 

Wing quarter-chord sucep angle, deg 

Table 2.1.1 

1.835 

r 

Ceonetric dihedral angle, deg 

Table 2.1.1 

7 

% 

Wing HEt coefficient 

Figure 1.1. 1.1 

variable 

Sun-aary : 

(Cy ) + (Cy ). - 2.32 x lCf 7 C 2 - 0.0007 per deg 

*5 7-0 ‘3 « 


< 














Symbol 

Description 

Reference 

Itafjnltude 

#; 

Wing-body poalclon parameter 

Figure 2.1.4 

0.326 

k i 

Wing-body interference factor 

Flyura 4. 1.2.1 

1.16 

2/3 

(C ) V 

‘ a f. 2 ' 3 
V 

C y of equivalent circular fuselage 

based on 2/3 power of fuselage volume 
(considered aqual to the negative of the 
potential flow part or (C T ) J tines 

a f- 2/3 

V 

the fuselage volume to 2/3 power 

Table 4.3.’ 
Reference 1 

0.002574 

S v 

Wing reference area, n 2 (ft 2 ) 

1 

Table 2.1.1 

14.40 (155.0) 

Sunnary: (Cy ), * 

T 3 * 

-0,002986 per deg 














Synbol 

Description 

Reference 

Magnitude 

a n 

Nuabtrof nacelle* 

* 

2 

S v 

Wing reference area, a 2 (ft 2 ) 

T ablc 2.1.1 

14 .40 (155.0) 

«*>■« 

Q 

Effective nacelle cross-section area, a.° (ft 2 ) 

Figure 2.1.4 

0.292 (3.14) 

in 

d a 

Effective nacelle fineness ratio 

Figure 2.1.4 

2.03 

(kj - k,) 

Nacelle reduced aass factor 

Figure 4. 1.3.1 

0.504 

Surxarv; (C, ) 
^3 n 

- -0.000476 per deg 















Symbol 

Description 

Rattrtnc* 

Magn-.cude 

s h 

Horizontal tail area, m 2 (ft 2 ) 

Table 2.1 

3.60 (3S.7) 

s 

V 

Vortical tall area, tu 2 (ft 2 ) 

Table 2.1 

1.85 (15.9) 

b v 

Vertical tall span, a (ft) 

Table 2.1 

1.52 (5.0) 

A 

V 

Vertical tail aspect ratio 

Table 2.1 

1.19 

(e A . 

Vertical tail chord at horizontal t?il, n (ft) 

Figure 2.1.S 

1.80 (5.92) 


Distance fron leading edge of vertical ta<l 
to a.c. of horizontal tail, in plane of 
horizontal tail, m (ft) 

Figaro 2,1.3 

1.45 (4.75) 

*c < c r> 
r h V 

Distance fron root chord of vertical till 
to root chord of horizontal ca.l. n (ft) 

Tlgure 2.1.5 

-0.051 (-0.167) 

<Ve 

Death of fuselage 3t quarter-rcoc cnord 
of vertical tail, a (ft) 

Figure 2.1.5 

0.661 (2.17) 

V £ > . 





Ratio of vertical tail aspect ratio in 

-Figure 4. 1,0,1 

1.5 

A v 

presence of fuselage to isolated vertical 
tail aspect ratio 

(a) 


A v (fh) 




X 1 

Ratio of vertical tail aspect ratio in 

Figure 4. 1.4.1 

‘ 1.25 

V*’ 

presence of fuselage and horizontal tail 
:o aspect ratio of tail in pretence of 
fuselage alone 

(b) 

1 

1 

. *h 

Relative tail size factor 

j Figure 4. 1.4. 2 

1.12 








































Table 4.1.5.1s L of the ATLIT Airplane 

• o 


Jynbol 

Description 

Reference 

Magnitude 

<S >v 
T a r-o 

Wing contribution without dihedral, 
pec dig 

Table 4. 1.1.1 

2.32xlO' ; 

V 

w 

(c ). 
U ‘ 

Increment due to dihedral, per deg 

Table 4. 1.1.1 

-0.0007 

«V' 

Fuselage contribution, par deg 

Table 4.1.:.l 

-0.00299 

(< V n 

Nacelle eoncrlbudon, per dee 

Table 4. 1.3.1 

-0.000476 

(C Y ^vfwfn) 

Vertical tail contribution with Inter- 
ference accounted for, per deg 

Table 4.1.4. l(c> 

-0.0056 

Sumnrv: (li, ) * -0.00977 per dee 

r prop 

J off 
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Propeller off lift curve of the wine 
for the ATI. IT airplane (Reference 1) 





























4.2 Yawing Moment Derivative, C n 

3 

The yawing moment due to sideslip derivative, C , of the complete 

n B 

airplane is found by considering the contributions of the following 
components s 

(1) Wing 

(2) Fuselage, including wing-fuselage interference 

(3) Nacelles 

(4) Vertical tail, including the interference and sidevash 
of the vertical tail. 


4.2.1 Wing Contribution to C 


The wing contribution to C is due to the increased induced drag 

n s 

on the leading wing caused by increased lift on the leading wing. 

Equation 4. 2. 1.1 from Reference 3 gives the incompressible (C n ) y 

3 



can A 


’V A v 


c M / 

co, A c/ .) ^ 


COS A , , “ T" 
c/4 


3 cos \ 


c/4 


o 


£ Sin A e/tV 

l n 


(4. 2. 1.1) 


where 

A w is the wing aspect ratio 

^ c / !t is the quarter chord sweep angle of the wing 
C y is the wing mean aerodynamic chord 

x is the distance (positive rearward) from the center of gravity 
to the wing aerodynamic center along the mean chord line. 




n 



To correct for the three-dimensional effects of compressibility, 

Equation 4. 2. 1.2 from Reference 3 gives the ving contribution to C . 

n„ 


J±) - * 4 e °‘ till \ / W + 4 Vi C0S Vc ~ 8 (%\ 

C lVm + 4 *°' •'«/*/ \ ^ + 4A U co* A c/4 - S cos“ A c/4 )( c ^) 


low 

speed (4. 2. 1-2) 


where 


B 2 - VI - M 2 cos-A c/4 


C L is the wing lift coefficient obtained from Figure 4. 1.1.1 


The summary calculations of the wing contributions to C are 

n 3 


shown in Table 4. 2. 1.1. For the ATLIT airplane 


( c n ) w =* 0.000156 C L 2 per deg 
3 w 


(4. 2. 1.3) 


4.2.2 Fuselage Contribution to C 

n. 

The fuselage plus wing-body interference contribution to C is 

n 3 

given «*.n Equation 4. 2. 2.1, taken troia Aeference 3* In general, wing 
sweep, wing taper, and Mach number have no effect. 


(S f )_ If 

(C n >£(„) ' -S "IT t 

D W W 


(4. 2. 2.1) 


where 


(S f ) s is the fuselage side area obtained fron Figure 2.1.4 


S is the wing area 
b is the wing span 


lc is the length of the fuselage from Figure 2.1.4 





.'i ' \C- ft 


t v,. 
k 


There are two ways to obtain the correlation factor K^. In 

Reference 4, is given for aidwing configurations and it is assumed 

•that (C^ does not vary with angle of attack. A nomograph for 

8 

obtaining midwing configurations is shown in Figure 4. 2. 2.1. 

In Reference 3 it was concluded on the basis of wind-tunnel data 
that angle of attack variations were significant for non aidwing con- 
figurations. Figure 4. 2. 2. 2 taken from Reference 3 13 used to extend 
the nomograph to obtain K^, based on angle of attack. This graph was 
constructed for an airplane with 2z w /(w^) w = 0.50 and is based on 
tail-off wind-tunnel data. For the ATLIT airplane ? - w /(v.) u “ 0.33; 
nevertheless. Figure 4. 2. 2. 2 was used to obtain in the absence of 
tail-off wind-tunnel data. If tunnel data is available, the procedure 
outlined in Reference 3 could be used to derive for other wing 
locations. 

Table 4. 2. 2.1 shows the summary calculation for the fuselage 
contribution to C . 

n s 

4.2.3 Nacelle Contribution to C 

n 3 

The nacelles contribution to C derives directly from the nacelle 

n 3 

side force acting through the nacelles moment arm. Equation 4. 2. 3.1 
from Reference 3 is used to obtain the nacelles contribution. 


■ x cos a + z sin a . 

(C n K - < C Y >n h > 

n_ n i * n b 

3 - 3 w 


(4. 2. 3.1) 


where 

(C„ ) is the nacelles contribution to side force due to side slic, 
Y 2 n 

from Section 4.1.3. 


0 / 



. V * * w* * / > J 

• 5? - 


is the longitudinal distar :e from the nacelle center of 
pressure to the airplane center or gravity, obtained from Figure 2.1.4 
is the vertical distance from the nacelle center of pressure 
to the airplane center of gravity, from Figure 2.1.4 
b y is the wing span, obtained froa Table 2.1. 

Table 4. 2. 3.1 summarizes the contribution of the nacelles to C 

n s 

4.2.4 Vertical-Tail Contribution to C 


The vertical tail contribution to C is obtained in a similar 

n 3 

manner to the nacelle contribution. The vertical tail C is the 

n s 

vertical tail C times the vertical tail moment arm and is obtained 
8 . 
from Equation 4. 2. 4.1 from Reference 3. 

/ cosa - z sina 

n^ v(wrh) Y^ v(vfn) b y . v 


. (C„ ) is the side force due to sideslip of the vertical 

Yg v(vrn) 

tail, from Section 4.1.4 

l is the distance from the center of gravity to the quarter 
chord of the mean aerodynamic chord of the vertical tail, parallel 
to the X-body axis, obtained from Figure 2.1.3. 

z v is the distance from the center of gravity to the quarter chord 
of the mean aerodynamic chord of the vertical tail, parallel to the 
2-body axis, obtained from Figure 2.1.3. 


Table 4. 2. 4.1 summarizes the vertical tail contribution to 





•/ 


4.2.5 C of Che ATLIT Airplane 
R 3 

The total yawing moment derivative of the complete airplane 
(propellers off) is 


(C n 3 } w + ^n^ffw) + ^n^n + 'SigVvfh) 


Table 4. 2. 5.1 summarizes the contribution of each component 
and the total C . 






N>. 









ruatltfto Contribution to 


Svrsbol 

Ascription 

Reference 

Magnitude 

<V. 

Fu-ioi^ia aid* area, n* tft‘> 

Tabic 4.3.1 
Reference l 

tsi.u 

3 U 

'din* reference are*. tar ift~) 

Table 2.1.1 

■ 

li.iii aJ5.ni 

/f 

lencto. ri (ft) 

Figure 2.1.4 

i.js: (r.40i 

b 

w 

Win* spin, 3 (ft) 

T.hle 2.1.1 

12.19 (40.05 

*« ■_ 

Verticil illstmce of wln< root chord 
below the* centerline o: the equivalent 
fuse line, n \:t) 

Figure 2.1.4 

0.104 vO.033) 


«ldth of eh** :*ua*l.i*e it the ;u.irtor 
chord ot the exposed root chorl, n <:t) 

ri<ure 2.I.* 

i.i9 o.?:> 

S *£ 

Revnold# mmo*r baaed on fuJ*la<e length 

- 

1.J7 x 10* 

* 3 *»• V h : 

ile one trie fuselico pirmeteri 

Figure 4.2.2. 1 

a» llated 

s 

Correct ton factor for fuiel.tc* contrL- 

ric-jre 

variable 

button 

' 

* . . .. .. 



r.*. i 


0 . Wl: 

d.WUJ 
VV,u * 


O.-Wl* 

v;o«u* 





















table 4. 2. 3,1: Nacelle Contribution to C 

"i 


Svr-iol 

Description 

Ref.rance 

Magnitude 

<S >„ 

Nacelle contribution to Cv 

Tibi. 4. 1.5.1 

-0.0004T6 


'i 



X, 

Distance rron c.g. to nacelle center 
of pressure, m (ft) 

Figure 2.1.4 

0. 966 0.17) 

r 

n 

Distance fron X-body axis to nacelle 
center of pressure, m (ft) 

Figure 2.1.4 

-0.1J2 (-.50) 

b 

V 

Wing span, n (ft' 

Tabla 2.1.1 

12.19 (40.3) 


^ I . Surctary: (C ) • -0. 0000377 cm i+ O.OOOOOoO ala a 

• n^ a 
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Table 4. 2. 5.1: C of the ATLIT Airplane 

** a 


a. 4eg 

V« 

Table 4.2. 1.1 

V' 

Table 4. 2. 2.1 

n j n 

Table 4. 2. 3.1 

(C } 

v(wt‘n) 

Table 4. 2. 4.1 

(C ) 
n prop 

i off 

Eq. 4. 2. 5.1 

_ | 

0 

-0.00013 

-0.00004 

0.00206 

0.20139 

-2 

0 

-0.30014 

-0.00004 

0.00208 

0.00130 


0.00001 

-0.00016 

-0.00004 

0.30210 

0.00191 

* 

0.00004 

-3.00017 

-0.00004 

0.00212 

0.00195 

4 

0.00007 

-0.00022 

-0.00004 

0.00213 

0.30194 











Axis of 
equivalt 
fuselagi 


| . 00 - 

o . 

I 


.00 







s 


a Degrees 

Figure 4.2.5-.1: Comparison of predicted C r with full scale 

: *3 

wind tunnel data (propellers removed, 

= 2.3 million) 
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if 


n 

I 
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n- 



4 • 3 Rolling Moment Derivative, C , 

'8 

The airplane rolling moment due to sideslip, C, , is composed 
of the following contributions: 

(1) Wing (without dihedral plus the effect of dihedral) 

(-) The effects of the fuselage on wing contributions 
(3) The vertical tail contribution. 

For general aviation airplanes without large dihedral angle 
on the horizontal tail, it can be assumed chat there is no significant 
contribution by the horizontal tail. 

4*3.1 Wing Contribution to C 

l 2 


From Reference 3 for subsonic speed in the linear lift, range, 
the wing contribution is given in Equation 4. 3. 1.1. 


\ 


‘V" ’ ( V"-„ * (c i>r * < c t >. 


“3 w r=o 1 -'*3 9 


(4.3.lil) 


Generally the effect of twist (C ) is not significant. 


The wirig contribution to C, . in the absence 

"3 

is given by Equation 4. 3.1.2 from Reference 3. 


of geometric dihedral 


V--0 = C L 

0 i =0 w 


C L / >1=0 
w> 


C L 

w 


where 


(4. 3. 1.2) 


C L is tfle win S lift coefficient from Figure i.l'.l.l 
w 

\\ 

C L / M=0 lb thS l0U SpGed ? ortion obtained from Figure 4. 3.1.1, 


M C , 

L 3 \ 

C L />! 1S the correction for compressible flow obtained f 


ron equation 


4.3. 1.3. 
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(4. 3. 1.3) 


where y* is the spanvise location of the centroid of the angle of 
attack load, as a function of senispan, obtained from Figure 4. 3. 1.2. 

The effect of unifrotn dihedral is accounted for by Equation 
4. 3. 1.4 from Reference 3. 


( % } r " r t-T* m=o \ ? er de s 


(4.3..1 


where 


T is the dihedral angle in degrees 


, (~ ~ i- 3 c h e incompressible effect of uniform dihedral obtained 
from Figure 4. 3. 1.3. 

is the compressibility correction from Figure 4. 3. 1.4 


The summary calculations of the wing contribution to C, are shown 

8 


in Table 4. 3. 1.1. For the ATLIT airplane 

(Cj ) w = -.000266 C L - .00154 
3 w 


(4. 3. 1.5) 


4.3.2 Effect of Fuselage on Wing Contribution to C, 

*3 

While the fuselage itself has negligible C . , it does provide 

4 a 

three wing-fuselage interference effects which are quite appreciable. 

The first effect is caused by the vertical position of the wing. 

The flow about the fuselage in sideslip induces vertical' flow components 




A -M 2 can A. 


iy . A + 4 f 2 j 2 * (_4_ y . A 2}l2 T7f j 

L\ cos A c /J u J I LV 03 'c/4 / a j 
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which cause a change in the effective wing angle of attack. Uithouc 

geometric dihedral, a high wing will produce a negative C. , a low wing 

B 

will produce Dositive C. , and a raid wing will generally show no change. 
Equation 4. 3. 2.1 from Reference 3 is used to obtain this first inter- 
ference effect. 


where 


^f(w) " 

*3 rw r-o 


1.2* / A~ 2 . 

w _w ( h + w) 

57.3 b " b 
w w 


(4.3. 2.1) 


is the vertical distance from the axis of the equivalent 
circular fuselage to the quarter chord of the exposed wing panels, 
positive down, obtained from Figure 2.1.4. 

h is the height of the fuselage at the wing location, from 
Figure 4. 3. 2.1 . - 

w is the width of the fuselage at the wing location, from 
Figure 4. 3. 2.1 


For wings with geometric dihedral, the position of Che wing with 
respect to the fuselage changes along the span. This is the second 
fuselage interference effect and is found by Equation 4.3.2 2 from 
Reference 3.; 

T per deg (4. 3. 2. 2) 

where 

(d r ) is the diameter of the equivalent circular fuselage at 

t w 

che wing, obtained from Figure 2.1.4. 

7 is the geometric dihedral, in degrees, obtained from Figure 
3.2.1. 


(C„ ) e ,^ s = -0.0005 >A 
*- 3 f (i ) w 
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V 


1 - 




o 

The. third fuselage interference effect is caused by the tendency 
of the fuselage to straighten the flow and therefore reduce the effective 
sideslip angle (see Sketch A). Reference 6 shows that this effect is 

I 

a function of forebcdy length, / f , wing span, and wing sweep. Ref- 
erence 4 indicates that the parameter A w /cosA c ^2 is also important. 

On the basis of Reference 6 this factor is considered zero for no 
sweep. For most general aviation' airplanes this is probably a reasonable 
assumption until systematic studies are made on this effect. 

Table 4. 3. 2.1 summarizes the fuselage interference contributions 

■ 7 ■ - 

to C. i of the ATLIT airplane. 

6 ' 



SKETCH A 

\ 

\ 


4.3.3 Vertical Tail Contribution to C - 

‘s ■ 

The vertical tail contributes to the airplane C. bv wav of the 
rolling moment produced by its side force due to sideslip. Equation 
4. 3. 3.1. from Reference 3 is used to findthe vertical tail contribution. 





























Table 4. 3.2.1: Fuselage Effect an Wing 


Synbol 

Description 

Reference 

Magnitude 

A u 

Wing aspect ratio 

Table 2.1.1 

10.32 

b 

u 

'./Ins span, a (ft) . 

Table 2.1.1 

12.19 (40.0) 

z 

w 

Vertical distance of wing root chord 
below centerline of thi equivalent 
circular fuselage, a (ft) 

Figure 2.1.4 

0.194 (0.63S) 

(d.) ■ h ■ w 

Dianeter of equivalent circular fuselage 
at wing, q (ft) 

Figure 2.1.4 

1.22 (4.0) 

r 

Wing dihedral angle, deg 

Table 2.1.1 

7 
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Vertical Tail Contribution to C 


Symbol 

Description 

Reference 

Magnitude 

^g’viufn) 

Vertical tail contribution to Cy 

3 

Tibia ft. 1.5.1 

-0.0056 

*v 

Distance perpendicular to X-body axis 
froa c.p. to quarter-chord o£ nean aero- 
dvnaaic chord of vertical tail, m (ft) 

Figure 2.1.5 

4.50 (14.75) 

fv 

Distance along X-body axis froa c.g. to 
quarter chard of mean aerodynamic chord 
cf vertical tail, n (ft) 

Figure 2.1.5 

1.22 (4.0) 

b w 

Wing span, a (ft) 

Table 2.1.1 

12.19 (40.0) 


Summary; (C t ) v ( u £ n ) ” 0.002065 sin a - n. 00056 cos a, per deg 


6 
























Figure 4. 3. 1.2: Spanwise location of centroid of angle 

of attack of loading for unswept wings 
(Reference 3) 





















4 


6 


a *v< Degrees 

Figure 4. 3. 4.1: Comparison of predicted C, with 

3 

full scale, wind tunnel data (propellers 
removed, N ' = 2.3 million) 




4.4 Snoiler Derivatives, C. , C 

*5 n S 

S3 

. In order to free Che entire trailing edge of the wing for high- 
lift devices, spoilers were used to provide roll control. The spoilers 
used are vented, gapped, upper-surface roll-control spoilers. For a 
discussion of the design decisions nade in developing the spoiler 
system, see Reference 5. Figure 4.4.1 shows the spoiler installation. 

4.4.1 Rolling Moment Due to Spoilers, C. 

u 

s 

The method used in obtaining the spoiler derivative C, is taken 

V 

s 

from Referenced and is based on a simplified lifting surface theory. 

This method is valid only for attached flow conditions. Reference 5 
outlines some of the characteristics of spoiler in the stall region. 

The rolling moment produced by deflecting one spoiler is given 
in Equation 4. 4. 1.1. 

c % 

Act' (4. 4. 1.1) 


where 

is the spoiler lift effectiveness in terms of the change in 
zero-lift angle of attack, obtained from Figure 4. 4. 1.1 

C! is Che rolling moment effectiveness of two full chord controls 
deflected antisymmetrically, obtained from Figure 4.4.1.? based on 
the effective inboard and outboard ends of the control. For partial 
span controls, this parameter is obtained as shown in Sketch A. 
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Sketch A 


The effective inboard and outboard spoiler locations are given 
in Equation 4. 4. 1.2 from Reference 4. 


n, = n, + i n, 
erf 1 1 


1 = n + A a 

°ef f 0 0 


(4. 4. 1.2) 


wnere 


r* i , t i o are the actual inboard and outboard ends of the spoiler 


given by n 


b/2 


u n^, A n o are the effective increments of the spoiler locations 
due to the spanwise flow of the spoiler wake, obtained from Equation 
4 . 4 . 1 . 3 . 


x; 

4(1 ~) cosA sin 6 

An. -;v [1 - (1 - X)r\ ] 


i A(1 + A) 


i cos(A TE + e)' 


(4.4.1. j) 


4(1 - -f) 


An 


o A(1 t .\) 


[1 - (5 - Dn 1 


cosA^sin 9 


o' cos(A T _, + 9) 



1 

1 
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where 

0 is determined from Figure 4.4.1. 3 as a function of spoiler 
sweepback given in Equation 4.4. 1.4. 

x s 

4(0.75 - (1 - “ ) ] , , 

tanA s - tanA c/A £ (i+T* (4. 4. 1.4) 

Using the above procedure, the spoiler rolling moment can be 

determined for various spoiler deflections and C. then found. 

<5 

s 

For the ATLIT spoiler deflections of 10, 20, 40, and 60 degrees 
were used. Table 4. 4. 1.1 summarizes the calculations for the ATLIT 
airplane. The predicted values are presented in Figure 4.4. ’.4 and 
compared with wind tunnel data. 





4.4.2 Yawing Moment Duo to Spoilers. C 

n 5 

s 

Tlu* yawing moment duo to spoilers comes from two effects. The 
first is due to a decrease in induced drag on one wing; the second is 
an Increase in profile drag on the wing with the deflected spoiler. 
The method used, taken from Reference 4, includes boch effects. 


For u.. swept wing at low angle of attack, the yawing moment is 
found directly from Figure 4. 4. 2.1. As in the ease of the rolling 
moment, should be computed for several spoiler deflection angles 








Mlt/i. 
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Figure 4.4.1: Spoiler installation (Reference 7) 



Table 4.4. 1.1: Spoiler Ef fectlveness of Che ATLIT Airplane 


! 

n 

l 


f 

I 

I 

I 

! 


Symbol 

Description 

Reference 

Magnitude 

ia ; 

Change In zero-lift angle of attack 
due to spoiler deflection, radians 

Figure 4.4. l.l 

variable 

S 

tolling oomenc effectiveness of two 
full chord spoilers, per rad 

Figure 4.4. 1.2 

0. 7507 

V n 0 

Inboard and oueboard nondlmenslonat 
spoiler locations, .j • 2y/b 

Figure 2.1 

0.421, 0.95 

n f • 
l eff °eff 

Effective spoiler locations 

Equation 4.4. 1.2 

0.410, 0.945 

3 

Spanwlse deflection of spoiler wake, 
deg 

Figure 4.4, 1,3 

-7 


Spoiler sweep angle. deg 

Equation 4. 4. 1.3 

-1.5 

4 , 

Spoiler deflection angle, deg 


variable 

\ 
■ c 

Distance from airfoil mean line to 
end of spoiler 

Figure 4.4.1 

variable 

X 

s 

c 

Distance froa airfoil leading edge 
to end of spoiler 

Figure 4.4.1 

variable 


4 J . deg 

X 

s 

c 

n ! 

da* 

s 

c t 

10 

0.771 

0.064 

0.056 

0.021 

20 

0.768 

0.0o6 

0.087 

0.035 

40 

0.755 

0.090 

0.108 

0.041 

60 

0.7J6 

0.110 

0. 123 

0.048 















p 


c 




Figure 4.4.1.2(c): Rolling moment effecnivenes 

(Reference 4) 
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6 , Spoiler deflection, deg 

igure 4.^.1. 4 Comparison or predicted spoiler effectiveness 
with full-scale wind tunnel data (Power off; 
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Table 4. 4. 2.1: Yawing Moment jue Co Spoilers on the ATLIT Airplane 


Symbol 


Erscrlptlon 


Reference 


Magnitude 


a 

c 


Spoiler span, a (ft) 

Wing caper ratio 

Distance from airfoil leading edge 
to end of spoiler 


Table 2.1 
Table 2.1 

Figure 4.4.1 


3.22 (10.53) 
0.5 

variable 


h 


Distance from airfoil mean line 
Co end of spoiler 


Figure 4.4.1 


variable 


S s , deg 

\ 

c 

h 

3 

C 

c 

n 

' 10 

0.771 

0.054 

0.0013 

20 

0.768 

0.066 

0.C027 

40 

0.755 

0.090 

0.0061 



60 


0.736 


0.110 


0.0127 


















t um my yyr ij 
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4.5 Rudder Control Derivatives. C„ , C . C, 

i » n > i * 

°R 5 R 5 R 

The rudder derivatives are found by assuming that the rudder 
is a simple flap. The side force is determined directly and the 
yawing and rolling moments are found by considering the side force 
acting through the appropriate moment arm. 

4.5.1 Side Force Due to Rudder 

The side force due to rudder deflection j.s found by Equation 

4. 5. 1.1 from Reference 3. 


< C L >v 
(C ) 

\ \ ' V’ 


(a 5 

r L 


(ct- ) r 

. °r C / 


*b 


(4. 5. 1.1) 


where 

(C^ ) v is the effective lift curve slope of the vertical tail, 
a 

obtained from Equation 4. 5. 1.2. This is the same as ~(C„ ) , ,, . 

Yj v(wfh) 

from Equation 4. 1.4. 5 with 3c/93 *» 0. 


a S 
v v 


(C i >» • K 1 (C L >„<fh) - s 

% a q w 

* m 


(4. 5. 1.2) 


(C. ) is the section lift curve slope of the vertical tail 
‘a 

from Reference 1,. Section 4.1 

. is the rudder span factor obtained from Figure 4. 5. 1.1 

\ , 

is the rudder chord factor obtained from Figure 4. 5. 1.2 


o'C, 

r l. 


. 4 r "/ J 


as a function of effective vertical tail aspect ratio and the ratio 

c c 

r r 

of rudder chord to vertical tail chord — . If — is not constant. , 

c c ’ 

v v 


an average value should be used. 


o 


is the section lift effectiveness of the rudder obtained from 
Equation 4. 5. 1.3 


\ ‘ “i ‘ C ( s > theory 


(4. 5.1. 3) 


where 


(C / 5 ) theory is the theoretical value of section lift effectiveness 


obtained fron Figure 4. 5. 1.3 

C / 

1 A 

is a correction factor obtained from Figure 4. 5. 1.4 


^ C l ^ theory 


as a function of c /c and 
r v 


r"', 


ri 


f a 

Tc 7T7 ’ wher - (C / J theory is obt3ined £ rom Equation 4. 5. 1.3 

/ a theory 1 a 7 

as a function of thickness ratio and trailing edge angle <f> T£ . 

(C / a ) theory “ 6.20 + 4.7(t/c) (1 + 0. 00375 <4^) per radio:. 


(4. 5. 1.3} 


k is a correction factor for reduced effectiveness at high 
deflection angles, obtained fron Figure 4. 5. 1.5. 

\ = •'l - M 2 

The summary calculations for C,„ of the ATL1T airplane are 

5 R 

shown in Table 4. 5. 1.1. Figure 4.5. 1.6 compares the calculations 
with full-scale wind-tunnel results for power off. 

4.5.2 Rolling Moment Due to Rudder, 

5 R 

The rolling moment due to rudder deflection is given by Equation 
4. 5. 2.1 from Reference 3. 


LOG 



/ 


r 

w 


O'i 






i v 


ot 


c, “ -c.. 

\ \ 


/' cosa + s' sina 
' v v 


(4. 5. 2.1) 


where 


/' is the distance parallel to the X-body axis from the center 
of gravity to the quarter chord of the nean aerodynamic chord of the 
portion of the vertical call spanned by the rudder. 

is the distance parallel to the Z-bodv axis from the center 
of gravity to the quarter chrod of the mean aerodynamic chord of the 
portion of the vertical tail spanned by the rudder. 

Note thac for full span rudders, and z^ are the same as / 
and z v defined in Section 4.2.4. 

Table 4. 5. 2.1 summarizes the ATI. IT airplane's C } . Figure 4. 5. 2.1 

V 

shows the comparison with full-scale wind-tunnel results. 

4.5.3 Va.'ing Moment Pm: to Rudder, C 

K 

The yawing moment due to rudder deflection is given by Equation 

4.5. 3.1 from Reference 3. ' 

/’ cosa - z 1 sina 
• v v 


C = -C 

n, \« 

■R K 


(4. 5. 3.1) 


where / * and ate defined in Section 4.5.2, 

* v v 

The surunarv caleulat ions for d of the ATI. IT aim lane are shown 

xv. 

a r 

in Table 4. 5. 3.1. In Figure 4. 5. 3.1 tae predicted values are compared 
with wind-tunnel results. 
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Table 4. 3.1.1: C„ far the ATLIT Airplane 

* 1 


Symbol 

Deecrlptlon 

Reference 

haqnltud* 

S v 


Vertical tall area 

Table 2.1.1 

1.35 (19.9) 

s 

w 


Wing Area 

Table 2.1.1 

14.40 (155.0) 

c l 


Effective lift curve slope of the 

Equation 4. 3. 1.2 

o 

IJ 

o 

VA 

\ 


vertical tall, per rad 



(c t K 


Section lift curve elope of the vertical 

Table 4. 2.4.1 

6.245 

‘a v 


call, per rad 

Reference 1 


s 


Rudder epan factor 

Figure 4. 5. 1.1 

1 


(a J 5 C, 






r L 
<*3 >C. 


Rudder chord factor 

Figure 4. 5.1. 2 

1.13 


r L 

V 




l 

V 


Vertical tail taper ratio 

Table 2.1.1 

0.459 

fr 

• 

Ratio of rudder chord to vertical tall 

Table 2.1 

• 3.20 

c 

• V 


chord 



B 1 


*1 - m : 


0.997 

c ,; 





1 J 

— 

Correction factor for lift effectiveness* 

Figure 4. 5. 1.4 

0.S3 

theory 





Theoretical section. lift effectiveness. 

Equation 4.5. 1.3 

3.60 

■4 


per rad 



K‘ 


Deflection anqle correction factor 

Figure 4.5. 1.5 

1.0 

°l 


Section lift ef fectlvencsa or rudder. 

Equation 4. 5. 1.3 

3.13 

r 


per rad 



Surmary : 

s * 

r 

1.00239 per. deg 







t 



n " b72 

Figure 4. 5. 1.1: Rudder span factor (Reference 3) 
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s 

I 


J 


Symbol 


lv 


Tabli 4. 5. 2.1: C, of the ATLIT Airplane 

r 


Description 


Change in side force with rudder deflec- 
tion, per deg 


Distance parallel to the X-body axis 
from the c.g. to the quarter-chord of 
the mean aerodynamic chord of the portion 
of the vertical tail spanned by the 
rudder, a (ft) 


Distance parallel to the Z-body axis from 
the c.g. to the quarter chord of the mean 
aerodynamic chord c* the portion of 
vertical tail spanned uy rhe rudder, tn (ft) 


Wing span, n (ft) 


Reference 


Table 4. 5. 1.1 


Figure 2.1.5 


Figure 2.1.5 


Table 2.1.1 


Magnitude 


0.002S9 


4.50 (14.75) 


-1.22 (-4.0) 


12.19 (40.0) 



S 

r • A 


r 


deg 

per deg 

j 

-4 

0.00036 


-2 

0.00031 • 

r /i*. 

0 

0.00029 

t; v - 

2 

0.00025 


4 

0.00021 


3 

0.00014 

' . -. . ' 

13.9 

-0.00001 

t • ' 


mil.' a 
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|* Symbol 




Tabic 4. 5. 3.1: 


Description 


of the ATDIT Airplane 


Change in side force with rudder deflec- Table 4.5.1. 1 
cion* per deg 


Figure 2.1.5 


Distance parallel to the X-body axis 
from the c.g. to the quarter-chord of 
the aean aerodynanic chord of the portion 
of the vertical tail spanned by the 
rudder, a (ft) 


Distance, parallel to the Z-body axis froa Figure 2.1.5 
the c.g. to the quarter chord of the nean 
aerodynanic chord of the portion of 
vertical tail spanned by the rudder, m (ft) 


Wing span, a (ft) 


Table 2.1.1 


Magnitude 


0.00239 
4.50 (14.75) 

-1.22 (-4.0) 
12.19 (40.0) 


a, deg 

\ 
per deg 

-4 

-0.00104 

-2 

-0.00105 

0 

-0.00107 

2 

-0.C01CS 

4 

-0.00108 

3 

-0.00110 

15.9 

-0.00110 











2 


Figure 4.5. 


a "v> Degrees 

.1: Comparison of predicted C with wind tunnel data 
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CHAPTER 5 

PREDICTION OF PROPELLER-OFF DYNAMIC DERIVATIVES 

In this chapter the propeller-off dynam<c (rate) derivatives 
are calculated. The methods of Reference 3 are used in all cases. 

The methods are based upon lifting-surface theory and are considered 
valid for attached flow conditions. It is assumed that the flow 
will remain attached up to near stall for conventional general 
aviation airplanes of moderate to high aspect ratio. All computed 
derivatives are in the stability axis system. 

5.1 Roll Damping Derivative, C , 

P 

The roll damping of the airplane is found by considering the 
contributions of the following components: 

(1) Wing-body 

(2) Horizontal tail 

(3) Vertical tail 

(4) Nacelles 

5.1.1 Wing-Body Contribution to C^ 

P 

For conventional general aviation airplanes with fuselage-width 

to wing-span ratios of 0.25 or less, the wing is the dominating factor 

and the fuselage may be ignored. For zero lift, C, of the wing nay 

> 

be found from Figure 5. 1.1.1. 

To account for dihedral and non-zero lift conditions, Equation 

5. 1.1.1 from Reference 3 is used. 



[v-L 


< C L 

a C T 


^ C i 

p 


- (C 3 ) 777 -^r— 77 T-f + (AC, ). „ (5. 1.1.1) 

0 s w c T 0 (C L ) W, =0 (C 1 ) r=o z p d * af? 


where 


(C. ) n is the zero lift C 0 obtained from Figure 5. 1.1.1 

X, W **U X, 

p c L p 

(C ) is the propeller off lift curve slope of the wing at the 
° °L 

lift coefficient being considered, obtained from Figure 4. 1.1.1. 

(C. ) n is the prooeller off lift curve slope of the wing St 

L W_ 

° C L 
zero lift 

(C i o )r 

* . is a correction factor which accounts for wing dihedral, 

; r»o 

p ■ 

obtained from Figure 5. 1.1. 2 

(AC, ) , is the drag-induced rolling moment caused by wing drag 
X. p drag 

in roll, obtained from Equation 5. 1.1. 2 from Reference 3. This incre- 
ment is small for high aspect ratios. 


«V >d„ 5 • - 1 

p 


L A + 2 cos (A !f ) , 

w , , „ -2 a u c/4 w 1 N 

•A .,) " *“ *«'* V + 4 «*<'./*>» ' •V* 

w c/4 w L J 

(5. 1.1.2) 


where 


(C_ ) is the wing zero lift drag coefficient, obtained from 
U w 
o 

Reference 1. 

The low speed wing-body contribution to is next modified to 

P 

account for non-zero Mach number according to Equation 5. 1.1. 3. 

r 1 K, + 4 cos ("_//.),. r 1 


r n A -4 COS v,.: r -1 

(C ) = — c - /4 'r 7 |(c ) 

l J vz A 3, + 4 cos (A I wf 

L P JM w 2_ c/4 w L p JM=0 
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where 

B 2 - /I " M^cos^A c/4 )~ 
w 

Table 5. 1.1.1 summarizes the calculations for the wing-body 
contribution to for the ATLIT airplane. 


5.1.2 Horizontal-Tail Contribution to 

P 

In general the horizontal tail contribution can be found in an 

analogous manner to the wing contribution. For normal tails it is 

of the order of one percent on the wing contribution. 

For zero dihedral tails in the linear lift region, Fquation- 

5. 1.2.1 (from Reference 3) gives the horizontal tail .referenced 

P 

to the wing area. It also accounts for wing interference effects. 


< c t >hf 




5. LV: h + 4 eo,( -'e/* , vJ l r \ 


(5. 1.2.1) 


where all quantities are referenced to the horizontal tail and defined 
in Section 5.1.1. 

In Table 5. 1.2.1 the horizontal tail calculations are summarized 
for the ATLIT airplane. 

5.1.3 Vertical Tail Contribution to C^ 

During rolling, the wing induces sidewash on the vertical tail 
due to unsymmetrical span loading. The equation used to determine 
the vertical tail contribution, frou Reference 3, is 

z v cosa + / v siting [ 2(c„cosa + / v sina) 

~ r Jr 

(5.1.3,; 


K = - 57 - 3(C i \ 



Hi 



where 


(C£ ) is Che vertical Call effective lift curve slope obtained 
a 

from Section 4.3.1. 

z v is the vertical distance parallel to Che Z-body axis from 
the c.g. to the vertical cail mean aerodynamic chord, obtained from 
Figure 2.1,3. 

/ is the horizontal distance (parallel to the X-body axis) 
frcm the center of gravity to the quarter chord of Che mean aerodynamic 
chord of the vertical cail, from Figure 2.1.3. 


3o 



is the rate of change of sidewash with wing-tip helix angle. 


There is no easy way to calculate 


3 


3a 



2 V 


This factor varies 


with angle of attack and is influenced by aspect ratio, sweepback, 

3m’ tail geometry. A value of .20 was used in Reference 3 and thus 
was also used here due to lack of better guidelines and the similarity 
of the airplane in Reference 3 with the ATLIT airplane. 

The summary calculations for the vertical tail contribution are 
given in Table 3. 1.3.1 for the ATLIT airplane. It can be seen Chat 
the effect of the tail is snail. 


5.1.4 Nacelle Contribution to C, 


The nacelle contribution to C, is caused bv me change in 

t 

P 

angle of attack induced by the roll rate.. Equation 5.1.4. 1 from 
Reference 3 gives this effect. 


(C. ) *= -114 . d (C, ) (: — ’>*• 

«. n i-» n o 

o i w 


(3. 1.4.1) 
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where 


y T is the distance parallel to the Y-axis from the thrust line 

to the center of gravity 

(C^ ) is the lift curve slope of the nacelles 
a 

The value of (C ) was found by caking the derivative of Equation 
l n 
a 

5. 1.4. 2 from Reference 1 (Equation 4.3.5). 

C L - 0.002031 a + 0.0000201 a 2 (5. 1.4. 2) 

n 

therefore 


(C. ) - 0.002031 + 0.0000402 a (5. 1.4. 3) 

L n 
a 

The nacelle contributions are summarised in Table 5. 1.4.1. 


5.1.5 C i of the ATLIT Airplane 
“■p 


The propeller-off C of the entire airplane is given by Equation 


5. 1.5.1. 


W ‘ <C » >«< * >hf ♦ «, ->v ♦ «h >n 
" - P P P 


(5. 1.5.1) 


P off p 


Table 5. 1.5.1 summarises the propeller-off C ; of. the ATLIT. 

‘‘P 

Figure 5. 1.5.1 presents C, .as a function of angle of attack. There 

P ’ 

is no wind tunnel data available for comparison. 
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Table 5*1. 1.1: UlnK~3ody Contribution t o C 


■ Symbol 

Description 

Reference 

'iagnltjdc 

M 

Mach nunbar 

- 

0.031’ 

"w 

f l-«l-coe\ a, 
C/4 


0.997 

A 

V 

Wine aspect ratio 

Table :.i.l 

10.32 

'c/4 

Wing quarter-chord sweep angle, deg. 

Table 2.1.1 

1.315 

X 

w 

Wing taper ratio 

Table 2.1.1 

0.5 

b 

V 

Wing span, m (ft) 

Table 2.1.1 

12.19 (40.0) 

d 

V 

Fuselage width at wing, n (ft) 

Figure 2.1.4 

1.22 (4.0) 

<>«. n 
P '• L *0 

Zero lift wing-body C. , per rad 

Figure 5. 1.1.1 

-0.325 

2* 

W 

Vertical distance from c.g. to quarter- 
chord of wing root chord, n (tt) 

. 

Figure 2.1.4 

0.25 (0.36) 

r 

Wing dihedral angle, dec 

Taole 2.1.1 

7.0 

cc i>r 

? 

Dihedral effect on 

P 

Figure 5. 1.1.2 

1.0 

- (C t ) r-o 
? 

‘ C L >v, 

1 C L-9 

Wing lift curve slope at zero lift, 
per rad 

Figure 4,1. 1.1 

5.09 


Wing life curve slope at ner rad 

Figure 4. 1.1.1 

5.09 

<S K 

Zero lift wing drag 

Table 1.12. ’..2 

0.0097 



Refer* ice 1 



(C, • -0. 32fc2 - 0.00)366 0, * per rad 

? w 













•’abol 


Description 


s h 

S v 

b h 


% 

q„ 


A.,- 

h 


‘ C < >h. 

P '-L.O 


< iC l >h 


p ‘ drag 


Sumnarv: (C. 1, . » 

i .if 


Horizontal tail area, nr (ft-) 

Wing area, nr (ft‘) 

Horizontal tail aspect ratio 
Horizontal tail span, n (ft) e 
Wing span, ni (ft) 

Dynamic pressure ratio at horizontal 
tail 

Horizontal tail quarter-chord sveep 
angle, dog 

2cro lift C, of the horizontal t3il 
*> 

Increment in C. dun to drag 

V 

p 

-0.0054:5 - 0.0001:0 c. 


per rad 


Reference 

Magnitude 

Table 2.1.1 

3.50 (33.7) 

Table 2.1.1 

14.40 (155.0) 

T bla 2.1.1 

4.75 

Table 2.1.1 

4.13 (13.56) 

Table 2.1.1 

12.19 (40.0) 

Table 4.9. 3. 3 
Reference 1 

1.0 

Table 2.1.1 

0 

Figure S.l.l.t 

-0.33 

Fqtiac.cn 5. 1.1. 2 

variable 
















i 
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Table 5. 1.3.1: Vertical Tall Contribution to C, 

A. 


Symbol 

Description 

Reference 

Magnitude 

<V>v 

Vertical tall effective life curve slope. 

Table 4. 3. 1.1 

0.00507 

1 3 ' 

per deg 



*v 

Vertical distance from Che c.g. to the 
vertical tali mean aerodynamic chord, 
a (ft) 

Figure 2.1.5 

-1.22 (-4.0) 

l v 

Horizontal distance from the c«g. Co 
the quarter-chord of the vertical tail 
mean aerodynamic chord, m (ft) 

Figure 2.1.5 

-4.50 (14.75) 

b « 

Wing span, a (ft) 

Table 2.1.1 

12.19 (40.0) 


Rate of change of sldevash with vtngtip 
helix angle 

- 

0.20 












Tibia 5.1. 4,1: Nacelle ContrlbuCi.cn to C 


Symbol 

Description 

Reference 

Magnitude 

< C L >n 
a 

Nacelle life curve slope, per deg 

Equaclon 5. 1.4.3 

variable 

y T 

Distance parallel to Y-axis iron thrust 
line to c.g. , a (ft) 

Figure 2.1.1 

1.90 (4.2) 

b 

u 

Wing span, m (ft) 

Table 2.1.1 

12.19 (40.0) 



(C 1 V per de * 

0l v (leg 

P 

Eq. 5.1.4 , 1 

-4 

-0.00515 

-2 

-0.00537 

0 

-0.00559 

* 

-0.00581 

4 

-0.00603 

15.9 

-0.00735 

















a, de: 


, deg 

* C t *v £ 

P 

Table 5. 1.1,1 | 

< C 4 

P 

Table 5.1.2. 1 

< C l >v 
P 

Table 5. 1.3.1 


(C lVop 

p cff 
Eq. 5.1.5 

-4 

-0.5273 

-0.0055 

-0.00186 

-0.00515 

-0.5398 

.2 

-0.5274 

-0.0055 

-0.00084 

-0.00537 

-0.5391 

0 

-0.5277 

-0.0055 

o.e 

-0.00559 

-0.5388 

2 

-0.5282 

-0.0055 

0,00965 

-0.00581 

-0.5289 

4 

-0.5290 

-0.0055 

0.00112, 

-0.00603 

-0.5394 

15.9 

-0.5359 

-0.0055 

-0.00030 

-0.00735 

-0.5491 
















eo.+i 







Figure 5. 1.1.1: Zero lift wing and horizontal tail , in 

P 

eluding fuselage effects (Reference 3) 
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a '■> Degrees 

Figure 5. 1.5.1: Predicted o£ the ATLIT Airplane 

P 

(propellers removed) 
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5.2 Yaw Damping Derivative, C 

r 

The yaw damping derivative is largely caused by the vertical tail 
and to some extent by the wing. The fuselage is of minor importance 
for conventional general aviation aircraft. 


5.2.1 Wing Contribution to 

A wing undergoing yaw rate will experience asymmetric lift distri- 
butions due to changes in velocity. These lift changes cause induced 
drag changes which produce a yawing moment. Equation 5. 2. 1.1 from 
Reference 3 gives the incompressible wing contribution. However, the 

profile drag term, (C_ ) , should be at the appropriate Mach number. 

u w 


(c ) 

n w 
r 


o 


’< C n >1 ’ 


'(1C, ) 2 ‘ 

r 


r 

c_ 2 

C L + 
w 

«~> >v 

L 


0 

w 


„ _ 


(5. 2. 1.1) 


where 


< C n >1 


— accounts for the C 
2 


w 

yawing of an isolated wing about 
Figure 5. 2. 1.1. 


due to induced forces resulting from 


its aerodynamic center, obtained from 


<“» ) 2 

r 

— 7 m — r— accounts for the unsvnr.etrical spanwise distribution or 

'w 

o 

profile drag, obtained from Figure 5. 2. 1.2. 

If the airplane c.g. does not coincide with the wing aerodynamic 

center, then the factor (C )./C, 2 must be corrected with the following 

n i l 
r v 


increment : 
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\( *‘°« A c/4 , A y \; t3 ° A c, 

C, 2 ‘ "| \ A « + 4 - 03A c/ 4 2coaA C /i/c w A » 

w - 

X 

, 4 c ° aA c /4 /;\2 ' V 1 

A “ + AC ° S '-‘=^ . UJ V C L 2 . 


A w + 4 co.A c/4 


where 


l C n 5 i"| 


(5.2.1. 2) 


■v M c/ ,-0 


x is the longitudinal distance from the c.g. to the wing aerodynamic 


center. 


< C n >1 
r 


is the zero sweep wing damping parameter obtained 


W JA c/4 -0 


from Figure 5. 2. 1.1 for A c ^»0. 

The general wing contribution equation is now 




'2 

+ tr \ I 

VW D 'w <C_ J 
0 D v 

o 


(5. 2. 1.3) 


* 1 

>. : j 

v : \ 

ir 4 


For the ATLIT the c.g. was assumed to be at the wing aerodynamic center. 

Table 5. 2. 1.1 summarizes the wing contribution ro C for the ATLIT. 

r 

5.2.2 Fuselage Contribution to C 


For conventional general aviation fuselages the fuselage contri- 
bution to C is small compared to the vertical tail. No analytical 
r 

design method is presented to calculate (C ).. 

r 

5.2.3 Vertical Tail Contribution to C 


*ne vertical tail is the orimarv contributor to C . Its efrect 

n 

r 

comes from a straight restoring moment due to the sideslip induced, 
by yaw rate and from the l3g of wing-fuselage sidewash, — — . * n ' 
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general there is no method to calculate the oscillating sidevash, 

and in fact Reference 3 suggests that it nay not be important. 

Assuming that the wing cancels the fuselage sldewajh, the vertical 

tal1 C n is 8lven by Ration 5. 2. 3.1. 
r 


(C n >v ‘ - 114 - S < c i >v(- 

r a \ 


' v cosa - Z'j sina 


where 

where (C^ ) , / y , and z are defined the same as in Section 5.1.3. 
a 

Table 5. 2. 3.1 shov;s the summary calculations for the vertical 

tail C p ' of the ATLIT airplane, 
r 

5.2.4 C r of the ATLIT Airplane 
r 

The propeller-off CL of the ATLIT airplane is given by Equation 
5. 2.4.1. 


^ C n ^prop “ ^ C n ^w + ^ C n + ' C n \ 


Table 5. 2. 4.1 summarizes the p-ooeller-of f C of the ATLIT 

' r • 

figure 5. 2. 4,1 shows as a function of angle of attack. Mo wind 

r 


tunnel data is available for comparison* 
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Reference 


HaEnttude 


Table 2.1.1 

10.32 

Table 2.1.1 

O.JO 

Table 2.1.1 

1.63$ 

- ' 

0 

Figure 5. 2.1.1 

-0.0143 

Flgur- 5. 2. 1.1 

- 

Figure S.2.1.2 

-0.239 

Figure 4. 1.1.1 

variable 

* 

Tabl« 4.12.1.2 
Reference 1 

o.oo?o: 












vertical Tail Concr l hut Ion to C 


Symbol 

Description 

Reference 

Magnitude 

(C l M v 

Vortical tail effective lift curve slope. 

Tabic 4.5. 1.1 

0.00507 

a 

per iieg 



r v 

Vertical distance from the c.g. to the 
vertical tail tean aerodynamic chord, 
m (ft) 

Figure 2,1.5 

-1.22 (-4.0) 

/v 

Horizontal distance from the c.g.’ to the 
quarter-chore of Che vertical tall ncan 
aerodynamic chord, a (ft) 

Figure 2.1.5 

4.50 (14.75) 

K 

w 

Wing span, a (ft) 

Table 2.1.1 

12.19 (40.0) 


_ > 


tC r. V p9r r * d 

t 

(e*. 5.:. 3. i) 


-O.J757 

-0.0774 

- 0 . 07^0 

-0.3304 

- 3.4416 

-0.0343 











































<n 


C: 


§ n 


1 

ia 


5 . 3 Roll Due to Yaw Rate Derivative, 

r 

Rolling moment due to yaw rate is mainly affected by the wing ' 
and vertical tail only for airplanes of conventional fuselage shape. 
Therefore only wing and vertical tail effects are calculated. 

5.3.1 Wing Contribution to C. 


The wing contribution to is produced by the spanwise lift 

r 

differencial produced by the yaw rate. For unswept wings, C„ is 

fc r 

linear with lift until C. . There is also an increment due to 

Sax 

geometric dihedral. For zero Mach number, the wing C,, is given 

r 

by Equation 5. 3. 1.1 from Reference 3. 


(C 


AC, 


M w 


c l Jr*o L w 
w . 


57.3 


(5. 3. 1.1) 


where 


_ , is the compressible change in C, due to wing lift 

l l / r=o r 

w / 

obtained from Equation 5. 3. 1.2. 


1 + 


L./r-o 


Vj ~ S 2 . ( Vi * : C03 A c/A 

; U W 3, 3 co a A e/4 ) \a u 0, - A cos i j C l r 

1 + * 2 003 A cM \ 

A W * 4 C03 ’-c/4 \ 3 / 


■M-0 


(5. 3. 1.2) 


where 


is obtained from Figure 5. 3. 1.1. 


L /r=M=0 
w ' 


B, = *'l - M-cosA.l . 

2 c/4 

The dihedral increment is obtained from Equation 5. 3.1.3. 
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AC„ 


_r _1_ 
“ 12 


ttA sin A 
w c/4 


.A +4 cos A . » f 
L w C/4J 


(5. 3. 1.3) 


The valid lift range of Figure 5. 3. 1.1 is dependent upon sweep angle. 

For conventional unswept general aviation airplanes, it is usable 

throughout the linear lift range. Also note that this procedure ignores 

the effect of center of gravity location. In general there does not 

appear to be enough data to account for this effect. 

Taale 3. 3. 1.1 summarizes the calculations of wing for the 

r 

ATLIT airplane. 

5.3.2 Vertical Tail Contributions to C, 


Excluding the effects of sidewash. 


3o 


rb 


, as in Section 5.2.3 


3 - 


w 

2v 


the aperiodic contribution of the vertical tail to can be obtained 

r 

from Equation 5. 3. 2.1 fron Reference 3. These calculations are sum- 
marized in Table 5. 3. 2.1 for the ATLIT airplane. 

/ z cos cl+ I sina\ // cosa - s si: 

, v ‘v W ‘ v v 


(C, ) =-114.6(C 


i 'v 
r 




b • ' 

w 


; ) 


(5. 3. 2.1) 


where (C^ ) , z , and are defined the sane as in Section 5.1.3. 


5.3.3 C of the ATLIT Airplane 
r 


The propeller-off C, of the ATLIT airplane is given by Equation 

r 


5. 3. 3.1. 


(C Z >prop * (C Z K + (C Z >v 
r off r r 


(5. 3. 3.1) 


Table 5. 3. 3.1 summarizes the propeller-off C, of the ATLIT. 

r 

Figure 5. 3. 3.1 shows as a function of angle of attack. 




o 


i 

1 

Table 5.3. 1.1: Wing Contribution to 


Syrabol 

Description 

Reference 

Magnitude 

M 

Mach number 

- 

0.018 

A v 

Wing aspect ratio 

Table 2. 1.1.1 

10.32 

X 

V 

Wing taper ratio 

Table 2.1.1 

0.50 

A c/4 

Wing quarter-chord sweep angle, deg 

Table 2.1.1 

1.835 

f 

Wing dihedral, deg 

Table 2.1.1 

7 

C L ■ 
u 

Wing lift coefficient 

Figure 4. 1.1.1 

variable 

/V) 

\ L u/'-M-0 

Wing contribution to C, for M»T*0 
r 

Figure 5. 3. 1.1 

0.240 

9 

Wing contribution to Cj for F«0 

Equation 5. 3. 1.2 

0.240 

■ 

Dihedral contribution to 

r 

Equation 5. 3.1. 3 

0.00074 


Suraiary: (Cj ) v - 0,*4 C L + 0.00074 per rad 


w 




Table 5.3. 2.1: Vertical Tail Contribution to C. 


Symbol 

Description 

Reference 

Magnitude 

< C l’>v 

a 

Vertical call effective lift curve slope, 
per deg 

Table 4. 5. 1.1 

0.00507 


Vertical distance from the c.g. to the 
vertical tail mean aerodynamic chord, 
n (ft) 

Figure 2.1.5 

-1.22 (-4.0) 

tv 

Horizontal distance froa the c.g. to the 
quarter-chord ot the vertical tall se&n 
aerodynanic chord, m (ft) 

Figure 2.1.5 

4.50 (14.75) 

b « 

Wing span, a (ft) 

Table 2.1.1 

12.19 (40.0} 


a, deg 

(Cj ) y , per rad 
r 

(Eq. 5. 3.2.1' 

-4 

0.0263 


0.0239 

0 

0.0214 

2 

0.0138 

4 

0.0161 

15.9 

-0.0011 




































5 . 4 Yaw Due 


CO Roll Rate Derivative . C 
- : — n 

p 

In this analysis only the effects of the wing and vertical tail 
are calculated, since they are the predominating contributors to C 

n * 
P 

5.4.1 Wing Contribution to C 

n 

P 

Wing rolling produces antisymmetrical lift loading which causes 
a yawing moment due to induced drag. Additionally roll-induced 
changes in angle of attack produce a change in viscous drag. Also 
dihedral effects are considered separate. Equation 5. 4 . 1.1 from 
Reference 3 gives the wing contribution. 



where 


' C n >1 



is 


r=o 


the zero dihedral contribution of the anti- 


symnetrical life and induced crag. This' term is 
obtained from Equation 5. A. 1.2 from Reference 3 which accounts for 
tip suction effects and center of gravity location. 


A * cos.'. 


:-«-o 


:/4 


■ w]. r ~ V ; 1 , , , - 

A v \ \ •* i \'~l i 


( 5 . 4 . 1 . 2 ) 


wnere 


(C n >1 

D 


is found f ror?. Figure 5 . i . 1 . 2 


A c/4^° 



To correct for compressibility effects, Equation 5. 4. 1.3 is used. 



. /•V 4eosA c/4 \/ A A*7 (A u B,*eo»r v/4 )tanJA c/4 ^ 

V.| 

k J 

r-o ' w 3 2 + '* cosA c/4/\ A u + v(A u + co»A c/4 )tin :, . e/4 / 

"1 

J* 

o 


r-M-o 


(5. 4. 1.3) 


where 

B_ » \1 - ii 2 cos 2 A 


2 ' ^ “c/4 

To account for dihedral, Equation 5. 4. 1.4 from Reference 3 is 


used. 




(5. 4. 1.4) 


w 

where (C^ ).,_q from Section 3.1.1 with compressibility effects 
P 

included. 

AC 

n 

The viscous drag parameter r— T^ i n Equation 5. 4.1.1 is found 


3a 


5C i 


from Figure 5. 4. 1.2. The second viscous drag term — — — is found by 

•la 1 

plotting viscous drag versus angle of attack, as shown in Figure 5.1. 1.3, 
and graphically determining the slopes. An analytical method is pre- 
sented in Appendix A. Reference 3 points out the importance of in- 
cluding the viscous drag contribution to C in obtaining good 

"p 

correlation vith wind- tunnel results. 

Table 5.4.1. 1 summarizes the wing contribution to C . 

n 

? 

5.4.2 Vertical Tail Contribution to C 


The vertical tail contribution is given by Equation 5.4.2. 1 from 
Reference 3, where all terms are as defined in Section 5.1.3., 
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(5.4.2. 1) 


I" 


O 


t 


-- \. 




* ..— *** * *. 




(C ) 
n v 
P 


-57.3CC. .) Mv c<330t + s v 

1 b 


sina 


2(2 v cosa+/ v sina) ^ 3a 


P b , 


w 

2V ■ 


Tabic 5. 4. 2.1 summarizes Che vertical tail contribution to C p 
for the ATLIT airplane. 


a 


r 


5.4.3 C of the ATLIT Airplane 
n 

P 

Considering the wing and vertical tail as the only significant 

contributors, is given in Equation 5. 4. 3.1. 

P 

(C „ Vop • >„ + \ 

p oi( ? p 

The calculations are summarized in Table 5.4. 3.1 and plotted as' 
a function of angle of attack in Figure 5. 4. 3.1. No wind tunnel data 
is available for comparison. 


• 


3 

3 

I 

5 1 o 

• 4 

.1 
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Symbol 


P«icript Ion 


Reference 


Magnitude 


M 

Mach cumber 

- 

o .am 

B, 

l-M*co*'\ ,, 

* 

0.997 


c/«* 



A 

w 

wing icpecc ratio 

TjjU 2.1.1 

10.32 

\ 

u 

King toper tula 

Table 2.1.1 

0.50 

'c;4 

Wi a* quarter-chord sweep .tn*Io # de* 

Table 2.1.1 

1.8)5 

b 

V 

Win* *?4P., a (ft) 

Table 2.1.1 

12.19 (40.0) 

K 

Win* aeon jorodynnmo chord, n (ft) 

Table 2.1.1 

1.225 (4.018) 

X 



0 

(c „ >1 

a.c. c*g. 



C. i 

2ero sweep d.iaptnq parameter 

Figure S.4,1,1 

-0.072 

W >0 




< C n > 




°L r . 

'oro dihedral dawptn* par.uetcr 1 

Kqujti on 5.4,1.) 

-0.0732 

u T m 'J 




AC 

n 

j> 

1 

• i 

Pihedral contrtlmf.loti to c 

| 

. 

r.ju.itlon 5.4,1. 4 i 

0.004: 

AO 

P 

! 


n 

Viscous 4ra< ;»ar.imt»t**r 

Figure 5.4.I.* 

2.059 

0 








H/: 

Viscous dr a* rite *>f change 

Ki euro 5 .4. 1.3 

| 

vtrlJbU* 












Tabid 5. 4. 2.1: Vertical Tall Contribution to C 

n 


Symbol 

Description 

Reference 

Magnitude 

<V>v 

a 

Vortical tall effective lift curve elope. 

Table 4. 5. 1.1 

0.00507 

par dog 



‘v 

Vortical distance from the' c.g. to the 
vertical call at an aerodynamic chord, 
m (ft) 

Figure 2.1.1 

-1 22 (-4.0) 

/v 

Horizontal llscance from tha c.g. to che 
. quarter-cho.-d of the mean aerodynamic 
chord of the vertical tall, a (ft) 

Figure 2.1.5 

4.50 (14.75) 

b v 

Wing span, n (ft) 

Table 2.1.1 

. 12.19 (40.0) 

.V 

Rat* of change of slJevash with wing tip 
helix anglu 


o 

IJ 

o 

— 


J, deg 

(C n ) v . per rad 
P 

(Eq. 5. 4. 2.1) 

—•4 

0.00198 

fc 

0.00272 

0 

0 

* 

-0.00280 

4 

-0.005*7 

15.9 

-0.023)3 
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Table 5.4. 3.1 


a. its 

V- 

(Tab la 5. 4. 1.1) 

-4 

0.00452 


-0.00884 

0 

-0.02219 

2 

-0.03355 

4 

-0.04690 

15.9 

-0.09534 


< C n>v 

P 

(Table 5.4.2. 1) 

< C „ > 

P Prop 
off 

0.00198 

0.00650 

0.00272 

-0.00612 

0 

-0.02219 

-0.00280 

-0.03635 

-0.00567 

-0.05257 

-0.02333 

-0.11367 






























i 







; ; :!wywM > . ' .:w *r 




CHAPTER 6 

PREDICTION OF POWER-ON AERODYNAMIC CHARACTERISTICS 

The effects of power especially from a propeller operation 
can have a significant effect on the aerodynamic characteristics of 
an airplane. Due to the highly interactive nature of the propeller 
slipstream, there are. usually a number of effects that must be con- 
sidered. Power effects are considered at three power settings: 
T'-O, T'-.0195, T'=.1970, where T'=Thrust/q S . 

6.1 Propeller Power Effects on Static Stability Characteristics 

For conventional twin engine airplanes the vertical tail is 
assumed to be outside the propeller slipstream. Therefore the 
propellers are assumed to produce no sidewash or dynamic pressure 
change at the vertical tail. 


6.1.1 Power Effects on C„ 


The contribution of power to Cy is given in Equation 6. 1.1.1. 

3 


<AC Y 3 ) power = (AC Y 3 ) M p + ^S^nUq) + (AC Yg ) n(a p ) 


( 6 . 1 . 1 . 1 ) 


where 

(AC^ ),j is the normal side force caused by the propellers 

3 ‘ P 

(AC ) — 

Ygn(Aq) is the change in nacelle C^ due to the increased 

3 

dynamic pressure 


n 


(*c v ) 


V 7 n(a ) is Che change in nacelle C v due to power i 
** ? * ** 


nduced 


sidewash at the nacelle. 
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The normal force contribution (AC y ) is given by Equation 

0 “ p 

6. 1.1. 2 from Reference 3. 


(iC Y m ~ # 
Y 0 N p 57 ‘ 


f / s /prop \ 


(6. 1.1. 2) 


where 


I ; 7*^3# 
& 


n is number of propellers 

f is the propeller inflow factor, obtained from Figure 6. 1.1.1 
Sp/prop is the disc area of each propeller 

^ C N )p is Che normal force parameter at T'»0 given by Equation 
a . c 

6. 1,1. 3 from Reference 3. 


( S>p- <S>0 

a a • 


K^-80.7 




( 6 . 1 . 1 . 3 ) 


is the empirical normal force factor of the propeller given by 


S ’ 262( />0.3R + 262 </>0.6R_ + - 35 <A.c 


0. 9R 
P P 


( 6 . 1 . 1 . 4 ) 


where (j^) is the ratio of blade width to propeller radius at the 
‘ P 

given propeller radius stations. 

[(C ) ] is the normal coefficient at K.,=S0.7 obtained 

ci . K n =80.7 ^ 

from Figure 6. 1.1. 2 as a function of blade ancle 3'. 

The change in C.^ of the nacelles due to increased dynamic pressure 
8 

is given bv Equation 6. 1.1.5 from Reference 3. 


''“ C Y, ) n(Ai) ” [ - (C V^n] 

2 d prop q 


(6 . 1 , 1 . 5 ) 


! I 
r \ 


l -X 

l ■} 
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a : jpiwj y r . 1 ”: * lyay » > ijwui^pt, JJji>;. L l w ^T^ T ggygy^^j y .W^ 


,o 


where 

[(AC^ ) ] is nacelle contribution with propellers off as 
8 off* 

determined in Section 4.1.3 

Aq n 

is the increase in dynamic pressure at the nacelle due to 


P; 


power given by 


iq n S w (T c /prop) 


( 6 . 1 . 1 . 6 ) 


irR ^ 
. P 


The. effect of propeller induced sidewash on the nacelle contri- 
bution is given by Equation 6. 1.1. 7 from Reference 3. 

' 3 o- V / Aq„ \ 


(AC Y A<o ) = “ t(AC Y Uprop V 3S 
8 ' P 8 off X 


1 + 


(6. 1.1. 7) 


where 


3a 


-rf is the change in prooeller induced sidewash with sideslip 
do 

given by 


30 

1# “ C 1 + C 2 <S >p ‘ 

a 


(6. 1.1. 8) 


where and C 0 are obtained from Figure 6. 1.1. 3. 

Table 6. 1.1.1 (a) - (c) summarises the effects of power on C„ 

3 • 

of the ATLIT -airplane. The power on C y is given bv 

3 

C Y - (C v ) + (AC y ) (6.1. 1.9) 

“» prop DOUBT 

S off. . 

Table 6.1.1. 1 (d) gives the power on C,. of the ATLIT. Figure 6. 1.1. 4 

'"3 

presents this data as a function of angle of attach and power setting, 

T 1 . Also the calculated values are coraoared ‘with wind-tunnel results. 

■ c ■ 



~ *y • \ ■' ^ / y ' /' ‘r?/Tp *.*‘~ >i- >i'j- ; j 

v ,./ .\s#%Q;y *>-rC«YX<b<*X ^ 

' • ' . '. 'V-‘. • - s'' ■ ■ "> -r. • ,' '.y\ /za -J-X.i 

?* l llt, , , WMW " ^ V1 "" **»■ -» i ■.» ■ iH 1 . *** + • * wj?**- ‘ * i*y g q? v ; w r » q i5_v w « ¥ ) r ' ~ «\^ 


O 


6.1.2 Power Effects on C 


r 


n 


«ffect of power on ls due co po „„. lndue| , d chaiiges ln ^ 


Its sate factors are assumed i„ t he po „,, t , nc „„ E ^ 

<«„ ) 

power V»_ ■*■ Ut 

6 “V'p 


: n 3 >po»cr ■ <\\ p + + («, ),„ , ■ (6. 1.2.1) 


me propeller uontal side force increment, ^ ^ ^ by 

Equation 6. 1.2. 2. 3 P 




x d cosa + sina 




(6. 1.2. 2) 


where 


“V»p 13 the InCre " ent ln C v “« to propeller penpal force. 

3 


from Section 6.1.1 


c.g. 


c.g. 


*p is the longitudinal distance fro. the propeller to the airplane 
t p is the vertical distance iron the thrust line to the airplane 


b w is the vin § span. 

Since the two nacelle increments. (AC ) and (ic j 

n j 3 n Uq) n 2 n ( a 0 > * 

• acc through the same moment arm. their efF^^r n u 

arm, tneir erfect can be considered together. 

where (6. 1.2. 3) 

(-iCy d )n <^> 3nd ( “% )n (o p ) 3re the ?cwor increments to C„ due to 
the nacelles, obtained from Section 6.1.1 5 

K„ is tho longitudinal distance trot, the nacelle center of pressure 


w0 che airplane c.g., from Figure 2 . 1.4 
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>*, im- 

'/ . •v'r' *- 

'•• - -<»-T ,!s- 


* T •*•'«*?* <ff *EtV 


z is thn vertical distance from the nacelle center of pressure 
n 

to the airplance c.g., from Figure 2.1.4. 

The total airplane C is 

3 


(C ) - (C ) „ + (AC ) 


n 'prop 
8 off 


n. power 

B 


( 6 . 1 . 2 . 4 ) 


Table 6. 1.2.1 presents the summary calculations for the ATLIT airplane. 
Figure 6. 1.2.1 gives total C as a function of angle of attack and 


power setting, comparing it with wind tunnel results. 


6.1.3 Power Effects on C, 

*3 

Two power effects on C. are considered. The first is caused by 

S 

the normal force produced by the propeller. 1 The second is the result 
of the lateral shift of the portion of the wing immersed in the propeller 
slipstream. Equation 6. 1.3.1 gives the total power increment. 


“Vp-" • <1C 'A + ( 4C AA> 


where 


(6. 1.3.1) 


(AC. ).. is the prooeller normal side force increment given bv 
Equation 6. 1.3. 2. 

/ -z cos tt + x sin 'x \ 

( \>3 ■ ( Vs r (6 - 1 - : >' 2 > 

2 p 3 p \ w / 

where all terms are the same as in Equation 6. 1.2. 2. 

When the airplane is yawed, the center of pressure of the wing 
shifts proportional to tan^S^s^.i. On one wing the shirr, will be inboard, 
while on the opposite wing the shift will be outboard. The net result 
is a rolling moment. Equation 6. 1.3. 3, from Reference 3, gives the . 
increment. 
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6.1.5 Spoiler Control 

For an airplane with ailerons, the eftect of pouer on roll 
control should be small. For spoilers It depends on the proximity 
of the spoilers to the propeller slipstream. For che ATLIT airplane, 
the portion of the spoiler Immersed in the slipstream Is quite small 
at low angle of attack, and none of the spoilor Is Immersed at 
higher angles of attack. (For data on the immersed span as a 
function of power and angle ot attack, see Table 5. 1.3.2 of Reference 1.) 
Tor this reason it was assumed chat the effect of power would be small. 

In addition there is no vind tunnel data that gives the isolated effect 
of power. The spoiler derivatives were therefore assumed not to vary 

with power setting. 






rc remen t ot C, Ou»? to Propeller Soraul For 


Symbol 

Deter ip t Ion 

Reference 

Ma<nitude 

a 

Number of propeller* 

- 

2 

% 

Propeller radius, m (ft) 

Tib Id :.i 

0.97 (J.17) 

% 

Width of propeller bled*, m (ft) 

* 

0.1:7 (0.417) 
it .3 R 

? 

0.157 (0.515) 
at .6 K 

P 



. 

0.108 tO.354) 
at .4 S 



• 

? 

Sp/pfop 

Propeller Il«c are., ,8 p“. <* J (ft*) 

- 

:.9) 01.57) 

S u (T'/ r rop) 





Propeller correlation parameter 

. 

v*' table 

8 r\ * 
? 




f 

Prone ller inflow parameter 

Figure 6. 1.1.1 

variable - 

K 

•n 

Propeller aide fore** parameter 

tsuaclon 6.1.1.-J 



Propeller blt.'.e Jo.; 


varial 1< 

<S ' ? 

Propeller normal force parameter 

Figure b.l.l.l 

variable 

-» ? v <J * 7 



'S 

t 

Propeller uorr.il force derivative 

Fquirioil *>• 1.1.1 

■ 

vi r table 














Table 6.1.1.UU)! Chance In Nacelle C„ Due to Incroa >ej Dvnanlc Pressure, UC V ) 



i 



Syrcbol 

Description 

Retcrence 

Magnitude 

<JC Y >n 

Propeller off nr. ells contribution to 

Table 4. 1.3.1 

-0.000476 

i prop 
off 

C Y , per Je* 
i 



fa 

Power laJuccd change in nacelle dynamic 

Equation 6.1. 1.6 

variable 

•’s 

preeaure 




1. UC '', , n(i5), per Je* 

U’q. 6.1. 1.0) ("q. 6. 1.1.5) 


-0 .. >00107 
- 1. 000230 



















lable 6.1.1.1(c): Change In "acelle C, Oue to Propeller Sldevaah, (iCy ) 


Symbol 

DeacrJptlon 

Reference 

Magnitude 

< iC Y >n 

Propeller off nacelle contribution to 

Table 4. 1.3.1 

-0.000476 

& prop 

C, , per deg 



off 




3<5 o 




£. 

3d 

Change In propeller sldcvaah with 
tldeallp 

Equation 6. 1.1. 8 

viriaole 

C l* C 2 

Propeller eldewaeh facr.ora 

Figure 6. 1.1.3 

variable 

4 ’n 





Povir Induced change tn nacelle dynanlc 

Equation 6. 1.1.6 

variable 

0- 

pressure 




-g 5. Per J * S 

(Eq. 6. 1.1. 6) (Eq. 6. 1.1. 6) (Eq. 6. 1.1. 7) 


0.000012 

0.000060 

0.000110 

















(Tjble 6.1.1.1(b)) 

<Vn<V 

8 P 

(Table 6.1.1.1(e)) 

(C Y * power 
on, 
per 

(Zq. 6.1.1 

0 

o.ocooi: 

-0.01046 

-0.000107 

0.000060 

-0.01062 


-o.ooo:3o 


0.000110 


-O.OICS) 
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vi Degrees 

Figure 6. 1.1. 4: Comparison of predicted power-on C.. with 

* 3 

wind tunnel data (?' = .1970) 
c 



r r 



Tabic 6. 1.2.1: Pouer-on C of the ATLIT Airplane 


SynboL 

Description 

Reference 

(iC T 

Propeller noraal side force derivative. 

Table 6. 1.1.1 

3 P 

per dee 


x„ 

Longitudinal distance frota the propeller 

Figure 2.1.4 

p 

to the c.g. , a (ft) 



Vertical distance from the propeller to 

Figure 2.1,4 

... p . 

the c.g., s (ft) 


b « 

Wing span, a (ft) 

Table 2.1.1 

U V n(i5> 

Change in nacelle C., due to Increased 

a 

dynamic pressure, per deg 

Table 6. 1.1.1 


Change In nacelle C.. due to prooeller 
3 

sidevash, per deg 

Table 6. 1,1.1 

x n • 

Longitudinal distance frcn nacelle center r 
of pressure to airplane c.g., a (ft) 

Figure 2.1.4 


Vertical distance free nacelle center of 
pressure to airplane c.g., a (ft) 

Figure 2.1.4 


Magnitude 


variable 
2.10 (6.59) 
-0.127 (-0.417) 

12.19 (40.0) 
variable 


0.97 0.17) 
-0.15 (-0.50) 


(i V > 7 


0.0915 0.1970 

■ 

0.0915 0.1970 

1 0 


^'"n.^prep 


0.001S9 -0.0001 -0.C001 -0.0001 

0.00190 -0.0001 -0.0001 -0.0001 

0.03191 .-3.0001 -0.0001 -0.0001 

0.00195 -0.0001 -0.0001 -0.0001 

0.00194 -0.0001 -0.0001 -0.0001 I 0 


(C ) 

n 

3 pover 


0.0915 0.1970 


-0.00001 0.00179 0.00179 0.00173 

-9.00001 0.00130 0.00130 0.00179 

-0.00001 0.00131 0.00131 0.00130 

-0.O0C01 0.00 l S3 0.00135 0.00134 

-9.00001 0.00134 Q.001’4 O.OOtSl 
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Table 6. 1.3.1: Power on C of Che A7LIT 
3 


Airplane 


Symbol 

Description 

P.eference 

< 4C y>h 

Bp 

Propeller normal side force derivative, 
per deg. 

Table 6. 1.1.1 

2 P 

Longitudinal .distance from the propeller 
to the c.ge, a (ft) 

Figure 2.1.4 

x p 

Vertical distance from the propeller to 
the c.g. , m (ft) 

Figure 2.1.4 

b w 

Ming span, n (ft) 

Table 2.1.1 

«P 

Distance along thrust line from the 
propeller to the wing quarter-chord, 
a (ft) 

Figure 2.1.4 

<ac t ) w(«5) 

Increase In lift due to power-induced • 
increased dynamic pressure 

Table 5. 1.3. 2 
Reference 1 

(4 Vv(« p ) 

Increment in lift due to power-inducea 
dovnvash 

Table 5. 1.3.2 
Reference 1 

if* 

36 , 

Power-induced sldewash derivative 

‘ 

Equation 6. 1.1. 8 


Magnitude 
variable 
2.10 (6.59) 
-0.127 (-0.417) 
12.19 (40.0) 

variable 

variable 


a. 'leg 

(C 2,>orop 
' off 

(AC l >N 
3 p 

<4 V(iq 

+ V 

^l^power 

on 



» i 

A c 

— » 
*c 

T c 



0 

0.C915 

0.1970 

0 

0.0915 

0.1970 

0 

0.0915 

0.1970 

“ «# ■ 

-1.00212 

\Q 

M3 

xo 

0- 

C.0001 

0.2002- 

-0.00212 

-0.002C2 

-0.00192 

.1 

-o.co:o<> 

xO 

v0 

vO 

0 

0.0002 

■ 0.C0C4 

-0.10209 

-O.OC1S9 

-0.00169 

0 

-0.00206 

-0.00001 

-0.00001 

-0.00001 

0 

0.0003 

0.0005 

-3.00207 

-0.00177 

-0.00157 


-0.C0203 

-0.00001 

-1.00001 

-0.00001 

0 

0.0003 

0.O007 

-3.00204 

-0.00174 

-0.20134 


-0.00200 

-0. 00002 

-0.00002 

-0.00002 

0 

0.0004 

0.003ft 

—3.00202 

— 1.00io2 

-0. CO 122 

5.* 

-7.00174 

-0.10004 

-0.00004 

-0.00004 

0 

0. 3004 

0.2003 

-0.00178 

-0.00133 

-0. CO 03 3 


j • 

i 
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6-2 Propeller Power Effects oh Dynamic Derivatives ' 

6.2.1 Power Effects on C„ 

£ 

P 

The power Increment on is caused by the propeller normal force, 

P 

the change in dynamic pressure over the immersed sections of the airplane 
and the change in downwash in the slipstream. The power-induced normal 
force contribution is given in Equation 6. 2. 1.1. 


( \2 

rj 

w/ 


( 6 . 2 . 1 . 1 ; 


where 


(C L ) N is the lift curve slope of the propeller normal force 
a p 

obtained from Figure 6. 2. 1.1 (from Table 5. 1.3.1 of Reference 1) 

y T is the lateral distance from the thrust line to the airplane 


center of gravity 

b is the wing span. 


The increment in wing due to power is given by Equation 6. 2. 1.2 


from Reference 3. 

(AC l p ) w(A5 + e p ) " (AC i p ) w(Aq) + (AC i r ) w( €p ) 


( 6 . 2 . 1 . 2 ) 


-11A.6 [(AC l + (AC l ) w(e } ] 

a ct p 


■v ,\ 2 


I 

\v 


where 


change in lift curve slope, due to the power- 
induced change in dynamic pressure obtained from Figure 6. 2. 1.2 (from 
Table 5. 1.3. 2 of Reference 1) 

^ AC L W) the chan 3 e - n curve slope due to the power- 

a P 

induced change in downwash obtained from Figure- 6.2.1. 3 (from Table 
5. 1.3. 2 of Reference 1). 
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The change in nacelle with power is given by Equation 6. 2. 1.3 

P 

from Reference 3. 

“V"^ ' (1C « P ) «W5) * (4C t p U p > 


<St >„ 


P prop 

off L 


q„. 


!Sl ' 

3a 


1 - 


3e 

u 

3a 


ft) 


(6. 2. 1.3) 


where 


(C^ ) is the propeller off nacelle contribution from 

p n prop 
off 

Section 5.1.4. , 

^ is the increment in dynamic pressure behind the propeller 

«- 

obtained from Table 5. 1.3. 2 of Reference 1. 


is the rate of change of propeller downwash with propeller 
“p 

angle of attack, obtained from Table 5. 1.3.2 of Reference 1. 

3t u 

is the propeller upwash gradient obtained from Table 5.1. 3.1 

of Reference 1. 

Total airplane is 
P 


" (C i ) prop + , Uc t ' ) w(q+e ) + (AC i ! ) n(Aq+« ) (6. 2. 1.4) 


off 


Tabic 6. 2. 1.1 summarizes the power effects on the ATLIT airplane. 


Figure 6. 2. 1.4 shows total as a function of angle of attack and 

P 

power setting. There is no wind-tunnel data available for comparison. 
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6.2.2 Power Effects on C 

. n 

r 

No analytical methods of computing the power increment to 

r 

were found. In Reference 3, the C variation with power was deduced 

r 

by comparing with wind-tunnel acta on a similar model. The same per- 
centage change in C with power for the model was assumed, 
r 

This procedure is highly empirical and cannot readily be extended 

to other aircraft. A comparison of the prediction of Reference 3 

with flight test data shows poor correlation. In view of this and 

due to the fact that no experimental ATLIT data is yet available, 

the power effect on C n is not predicted for the ATLIT. 

r 

6.2.3 Power Effects on 

r 

Reference 3 shows chat the effect of power on is due mainly 

r 

to wing contributions and that even this effect is negligible (of the 

order of 1™). For this reason the power effects of are not 

r 

computed for the ATLIT and are assumed zero. ' 

6.2.4 Power effects on C 

n 


Due to a lack of design data and for the same reasons as those 
outlined above, the power effect on C for the ATLIT airplane is 

"p 


assumed zero. 
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i j 

Table 6. 2. 1.1: ?over-on C t of the ATL1T Airplane 


• . -J 

1 

[ " 


p 




* 

Symbol 

Description 

Reference 

Magnitude 

> '.U ! 

\v 

^ C L *N 
a p 

Life carvo slope of Che propeller 

Figure 6. 2. 1.1 

variable 

. • *( 

[ ‘ 

normal force 



H' i :.t 

: r 

*T 

Lateral distance from thrust line to 
airplane c.g., m (ft) 

Figure 2.1.1 

1.39 (6.21) 

r :‘:i 

\ 

, ■ • 


Wing span, a (ft) 

Table 2.1.1 

12.19 (40.0) 

'.’j 


(iC L J ) w(45) 

Change In wing lift curve slope iue to 

Figure 6. 2. 1.2 

variable 

;^.Vj 


power-induced change in dynamic pressure | 



• V'l 

• * 

(4C L ’-(cl 
a p 

Change in wing lift curve slope due to 

Figure 6. 2. 1.3 

variable - 

’ 


propeller downwash 






Change in dynamic pressure behind the 

Table 5. 1.3. 2 

variable 




propeller 

Reference 1 



• . . 

(C 1 >n 

Propeller off of the nacelles 1 

Table 5. 1.4.1 

variable 


i 

P prop 

P | 



- ’!■ 

'(f 

i . 

% 

i 

off 

3e o 




'* ~ V '» 

»•] 


3 a 

P 

■ 

3c 

u 

3a 

Rate of change of propeller downwash 
with propeller anglo of attack 

. i 

Propeller upwash gradient 

Table 5. 1.3. 2 
Reference 1 

Table 5. 1.3.1 
Reference 1 

-.155 




<C l Vop 

P off 

(Table 5. 1.5.1) 


(iC t \ 
P P 


(Eq. 6. 2.1.1) 


(iC l + c p ) 


(Eq. 6. 2. 1.2) 


-0.5398 

-0.3391 

-0*5338 

-0.5289 

-0.5394 

-0.5491 


-0.00099 

-0.00110 

-0.00121 

-0.00099 

-0.00110 

. -0.00121 

-0.00099 

-0.00110 

-0.00121 

-0.00099 

-0.00110 

-0.00121 

-0.00099 

-0.00110 

-0.00121 

-0.00099 

-0.00110 

-0.00121 


0 

.0915 

.1970 

0.00083 

-0.0153 

-0.0357. 

0.00083 

-0.0152 

-0.034C 

0.00083 

-0.0152 

-0.0284 

0.00083 

-0.0152 

-0.0284 

0.00083 

-0.0105 

-0.0284 

-0.00276 

J 

0.0028 


(iC l p ) n(45 + c p ) 
(Eq. 6. 2.1.3) 


(C, ), per rad 


(Eq. 6. 2. 1.4) 


0 

.0915 

.1970 

0 

.0915 

.0970 

0.00011 

-0.00059 

-0.00146 

-0.5399 

-0.5567 

-0.5792 

0.00012 

-0.00062 

-0.00153 

-0.5391 

-0.5560 ' 

-0.5758 

0.00012 

-0.00064 

-0.00159 

-0.5388 

-0.5557 

-0.5700 

0.00013 

-0.00067 

-0.00165 

-0.5289 

-0.5459 

-0.5602 

0.00013 

-0.C0069 

-0.00171 

-0.5394 

-0.5517 

-0.5707 

0.00016 

-0.00085 

-0.00209 

-0.5527 

-0.5511 

-0.5496 
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a, deg 

Figure 6:2. 1.1: Increment in lift due to propeller 

normal force (F.eference 1) 















Figure 6. 2. 1.3: Increment in lift due to propeller 

dcwnwash (Reference 1) 
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a 'v Degrees 

Figure 6. 2. 1.4: Effect of power on predicted of the 

P 


ATI.lT Airplane 






CHAPTER 7 

CONCLUSIONS AND RECOMMENDATIONS 



This report presents an analytical method for predicting the 
lateral-directional aerodynamic characteristics of light twin- 
engine airplanes. This method is applied to the Advanced Technology 
Light Twin (ATLIT) airplane, and the predictions are then compared 
to full-scale wind tunnel data. 

Based upon this comparison, the following observations and 
conclusions are presented: . 

1. The predicted values of C y arc rlwn the wind 

B 

tunnel values; however, the overall value is fairly 

reasonable, especially for preliminary design. It is 

necessary to consider more completely- the effect of 

angle of attack. The two areas that <-re thought to 

need the most study are the fuselage C Y variation 

8 

with angle of attack and the fuselage-vertical tail 
interference variation with angle of attack. 

2. For the power-off case at least, the average predicted 

values of C n show fairly good agreement with wind 

tunnel data. The tunnel tests show large variations 

of C with angle of attack,' especially with power on. 
n B 

The contribution of the fuselage Is strongly influenced 
by the wing location, and only partial methods are 
available to account for it. A more detailed study 
needs to be conducted to determine the effect of wing 
location and angle or attack." 
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3. The predicted C. of the ATLIT is slightly higher than 
wind tunnel results but still compares quite well. 

4. The spoiler performance, both rolling and yawing moments, 
is predicted somewhat higher than found in the wind 
tunnel data. However, the agreement is still quite 
good considering the highl” general nature of the 
methods used. A method should be developed which is 
tailored to general aviation airplanes. 

5. The prediction of the rudder derivatives was dis- 
appointing. But again it appears that the angle 

of attack factors have not been accounted for. The • 
angle of attack variation of the wing and body wake 
needs to be included in the prediction method. 

6. The dynamic derivatives are also presented in this 
report. When experimental data becomes available, 
it should be correlated with these predictions to 
evaluate the prediction method. • 

i. There needs to be more power-on wind tunnel data 
in order to better correlate the effects of power. 

Only limited power data was available, and it was 
at a relatively low thrust coefficient. 

8. The data in this report and in Reference 1 should be 

correlated with flight test data as it becomes available. 
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APPENDIX A 

A SIMPLE METHOD FOR COMPUTING 

MING VISCOUS DRAG, C p' 

o 
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APPENDIX A 


A SIMPLE METHOD FOR COMPUTING 

WING VISCOUS DRAG, C„' 

0 

o 

The wing drag which varies with angle of attack is comprised 

o’ two components. First is the drag due to direct lift production 

(vortex drag). The second is a type of profile drag caused by 

the buildup of the upper surface boundary layer as angle of attack 

increases (viscous drag). Together these two factors comprise 

the induced drag. Since the viscous drag term is needed for the 

prediction of , the following procedure from Reference 6 is 
P 

given. 

Equation B.l gives the viscous drag of a wing. 


C D ' - k A (A.l) 

o 

where 


k is an empirical sweep correction factor obtained from 
Figure A.l as a function of sweep and the parameter 
J siven in Equation A. 2 

A is the empirical viscous drag increment factor from 

Figure A. 2, also a function of J 

# 

and where 


\ 


°.3(ci + 1) — cos , ^ l 


(c. + l)(c,'+ 1) - 


(c., +■ l)tnn,j e 


(A. 2) 


where 


c^ and c., are taper ratio constants obtaine'd from Figure A. 3. 
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